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ABSTRACT 


Function-theoretic  techniques  which  have  been  successful  in 
solving  electromagnetic  scattering  problems  involving  wedges  satisfying 
an  impedance  boundary  condition  are  used  to  formulate  the  scattering  by 
a  wedge  composed  of  two  resistive  sheets.  The  goal  is  to  investigate 
the  possibility  of  generalizing  the  techniques  to  the  two  region  problem, 
characterized  by  nonzero  fields  in  both  the  Interior  and  exterior  of  the 
wedge. 

The  particular  techniques  considered  are  the  method  of  Maliuzhlnets' 
(1950  et  seq.)  and  the  Kontorovich-Lebedev  transform  (1939).  Both 
techniques  Involve  Integral  representations  (transformations)  for  the 
unknown  fields,  and  generation  of  functional  equations  for  the  transformed 
unknowns.  It  is  shown  that  for  the  resistive  wedge,  the  presence  of  an 
interior  field  results  in  functional  equations  of  considerably  greater 
complexity  than  those  encountered  in  similar  single  region  problems,  and 
which  cannot  be  solved  using  standard  techniques. 

A  novel  procedure  is  developed  which  replaces  the  functional 
equations  with  Integral  equations  of  the  Fredholm  type,  which  are  well 
understood  in  the  literature.  Though  exact  solutions  are  not  found, 
an  approximate  solution  Is  obtained  via  the  method  of  successive 
approximations  from  linear  operator  theory.  The  approximation  consists 
of  a  power  series  In  terms  of  the  resistivity,  and  bounds  for  its 
region  of  convergence  are  given. 
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CHAPTER  I.  INTRODUCTION 


1.1  Previous  Analyses  of  the  Scattering  from  Wedges:  A  Review 

The  scattering  of  electromagnetic  waves  by  wedge-shaped  regions 
has  been  investigated  by  several  authors,  although  few  exact  solutions 
exist  [1-28],  For  the  case  of  a  perfectly  conducting  wedge,  the  first 
solution  is  attributed  to  Somraerfeld  (see  Carslaw,  [1]),  and  to  Mac¬ 
Donald  [2].  These  works  addressed  the  two-dimensional  scalar  problem, 
and  were  based  on  an  extension  of  the  method  of  images  described  by 
Carlslaw  [3].  A  special  case,  in  which  the  included  angle  of  the 
wedge  is  zero,  was  elegantly  solved  by  Sommerfeld  [4],  This  solution 
for  the  half-plane  has  served  as  a  model  for  investigating  the  effects 
of  edge  diffraction.  Clemmow  [5]  developed  a  technique  for  general¬ 
izing  these  results  to  the  case  of  obliquely  incident  plane  wave 
excitation.  Modal  expansions  for  the  solution  can  be  found  in  [6], 
for  example,  while  an  expression  for  the  general  vector  problem  in 
the  form  of  a  dyadic  Green's  function  is  given  by  Tai  [7]. 

In  the  event  that  the  wedge  is  not  perfectly  conducting,  an 
exact  solution  cannot  be  found.  However,  for  the  case  in  which  the 
conductivity  Is  large  but  finite,  Leontovich  [8]  has  developed  an 
approximate  boundary  condition,  known  as  an  impedance  boundary 
condition,  which  allows  several  otherwise  intractable  problems  to  be 
solved. 
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The  major  attribute  of  this  approximation  is  that  it  replaces 
the  fields  within  the  conducting  body  with  approximate  surface  currents 
via  the  boundary  condition.  With  the  interior  field  thus  accounted 
for,  the  number  of  unknowns  is  halved,  and  only  the  exterior  field 
need  be  determined.  In  particular,  Senior  [9]  found  a  solution  for 
the  Impedance  half-plane  under  normal  plane  wave  incidence,  and  later 
extended  it  to  include  oblique  incidence  [10].  This  extension  was 
related  to  Cleramow's  technique  by  Williams  [11].  The  results  are 
obtained  through  application  of  the  Wiener-Hopf  method  (see  [12])  to 
the  unknown  currents  excited  on  the  half-plane.  For  a  wedge  of 
arbitrary  angle  and  differing  face  impedances,  a  general  solution  to 
the  scalar  problem  with  plane  wave  incidence  was  developed  by 
Maliuzhlnets  in  his  doctoral  thesis  and  subsequent  works  [13-15]. 

His  method  is  a  further  generalization  of  the  method  of  Images, 
and  involves  an  integral  representation  for  the  field  along  a 
Sommerfeld  contour.  Maliuzhinets  replaces  the  integral  equation  for 
the  unknown  by  an  equivalent  functional  difference  equation.  Senior 
[16]  and  Williams  [17],  following  similar  methods,  also  arrived  at 
the  solution  for  the  case  in  which  the  face  Impedances  are  the  same. 

Variations  on  the  half-plane  problem  and  its  solution  are  given 
by  Rawlins  [18],  Hurd  [19],  and  Bucci  and  Franceschetti  [20],  which 
are  essentially  founded  on  either  the  Wiener-Hopf  or  Maliuzhinets' 
method. 

Another  function-theoretic  technique  for  solving  scattering 


problems  In  wedge-shaped  regions  was  presented  by  Kontorovlch  and 
Lebedev  [21]  in  the  form  of  an  integral  transform  bearing  their 
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names.  A  solution  for  the  perfectly  conducting  half-plane  was 
given  as  an  example.  In  a  recent  paper  [22],  Jones  derives 
sufficient  conditions  for  the  existence  of  the  transform  and  a 
generalized  inverse.  In  a  major  extension  of  their  work,  Lebedev 
and  Skal'skaya  [23]  applied  the  transform  to  the  impedance  wedge  and 
developed  a  closed  form  solution  for  a  special  set  of  wedge  angles. 
Their  method  resulted  in  a  functional  difference  equation,  paralleling 
Mai iuzhi nets'  technique,  but  the  equation  is  of  second  order  instead 
of  first.  In  fact,  Maliuzhinets  [24]  described  the  relationship 
between  his  method  and  the  Kontorovich-Lebedev  (K-L)  transform  via 
a  Fourier  transformation. 

In  addition  to  the  impedance  boundary  condition,  of  which 
infinite  conductivity  and  its  dual  are  special  cases,  another 
approximation  exists  which  accounts  for  the  material  properties  of  a 
scattering  body  via  an  equivalent  boundary  condition.  In  this  case, 
the  body  is  assumed  to  consist  of  a  thin  dielectric  shell  of  a 
particular  shape,  and  is  approximated  by  an  infinitely  thin,  partially 
transparent  layer.  Because  it  is  most  applicable  to  lossy  dielectrics, 
this  approximation  is  referred  to  as  a  resistive  boundary  condition, 
although  due  to  the  transparency  of  the  sheet,  and  hence  the  existence 
of  both  Interior  and  exterior  fields,  it  is  more  exactly  a  transition 
condition,  analogous  to  those  applied  at  an  interface  between  two 
dielectrics.  Derivations  of  this  condition  can  be  found  in  [25,26]. 
Very  few  exact  solutions  for  resistive  bodies  exist.  Senior  [27]  and 
Anderson  [28]  give  solutions  for  a  resistive  half-plane  using  a  Wiener- 
Hopf  analysis.  An  excellent  review  of  the  impedance  and  resistive 
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boundary  conditions  for  the  half-plane  is  presented  by  Senior  in 
[29],  and  a  condensed  version  comparing  the  two  is  contained  in 
Section  1.3  of  this  work. 

Finally,  it  should  be  noted  that  the  electromagnetic  problems 
discussed  above  have  analogs  in  other  disciplines  involving  wave 
phenomena,  especially  in  acoustics.  The  acoustic  equivalents  of 
the  impedance  and  resistive  boundary  conditions  in  the  context  of 
Babinet's  principle  are  pursued  by  Senior  in  [30]. 

1 .2  Outline  of  the  Scattering  Problem 

In  order  to  further  understand  the  nature  of  the  scattering 
from  partially  penetrable  objects,  in  particular  those  composed  of 
resistive  materials,  an  investigation  into  the  scattering  of  an 
arbitrarily  polarized  electromagnetic  plane  wave  normalTy  incident 
on  a  wedge  with  resistive  faces  is  presented  here.  The  resistivity 
of  the  wedge  is  assumed  to  be  a  complex  scalar  constant,  independent 
of  position  on  the  wedge.  The  unknown  fields  exist  in  two  regions, 
the  exterior  of  the  wedge  (Region  1),  from  which  the  plane  wave  is 
incident,  and  the  interior  of  the  wedge  (Region  2),  as  shown  in 
Fig.  1.1.  The  primary  distinction  between  this  problem  and  those 
discussed  previously  is  the  existence  of  a  nonzero  field  In  the 
interior  region.  This  field  severely  complicates  the  task  of  finding 
a  solution.  In  some  respects,  the  problem  closely  resembles  the 
scattering  by  a  dielectric  wedge,  where  an  Interior  field  Is  also 
present.  Indeed,  exact  solutions  for  dielectric  bodies  are  very  few, 
and  those  which  have  been  developed  for  special  cases  of  the 


gure 
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wedge  [31-32],  are  either  extremely  complicated  or  somewhat  In 
doubt  [33]. 

Following  a  brief  derivation  and  comparison  of  the  Impedance 
and  resistive  boundary  conditions  in  Section  1.3,  a  formulation  of 
the  scattering  problem  to  be  considered  is  presented  in  Chapter  II. 

By  appropriately  decomposing  the  incident  field  into  two  components, 
with  either  the  E-  or  H- field  vector  parallel  to  the  edge  of  the 
wedge,  the  solution  can  be  constructed  from  a  pair  of  scalar  two- 
dimensional  fields.  The  problem  is  further  simplified  by  decomposing 
the  scalar  fields  into  components  which  are  either  symmetric  or 
antisymmetric  about  a  plane  bisecting  the  wedge.  Thus  a  total  of  four 
quantities  are  required  to  specify  the  general  solution. 

The  approach  taken  in  this  work  is  to  examine  the  feasibility 
of  extending  the  function- theoretic  techniques  used  in  solving  the 
impedance  wedge  scattering  problem  to  the  resistive  wedge.  While  the 
Wiener-Hopf  method  has  successfully  been  applied  to  the  half-plane  and 
other  geometries  (see  for  example,  Carlson  and  Heins  [34,35])  it  is  not 
appropriate  for  wedges  of  arbitrary  angles,  with  the  possible  exception 
of  the  right-angle  wedge.  Senior  discusses  this  shortcoming,  as  well 
as  certain  analogies  between  Malluzhinets'  method  and  the  Wiener-Hopf 
technique  in  [16,29].  Therefore,  the  emphasis  here  will  be  on  the 
methods  of  Maliuzhinets  and  the  Kontorovich-Lebedev  transform.  It  is 
shown  that  both  methods  lead  to  difference  equations,  as  was  the  case 
for  the  impedance  wedge.  However,  when  applying  the  method  of 
Maliuzhinets,  the  resulting  difference  equation  for  the  resistive  wedge 
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is  of  third  order,  compared  to  the  first  order  equation  obtained 
for  an  impedance  wedge.  In  general,  a  straightforward  technique 
for  solving  this  third  order  equation  is  not  available  [36],  A 
derivation  of  this  equation  and  its  consequences  is  given  in 
Chapter  III. 

In  Chapter  IV,  the  relationship  between  the  representation 
of  Maliuzhinets  and  Kontorovich  and  Lebedev  is  presented  [24]. 

The  latter  is  applied  to  the  resistive  wedge  formulation,  and  again 
a  set  of  difference  equations  for  the  unknowns  is  generated. 

As  before,  the  equations  do  not  yield  to  the  same  methods  of 
solution  available  for  the  impedance  wedge.  However,  a  novel 
technique  for  converting  the  difference  equations  to  Fredholm 
integral  equations  of  the  second  kind  is  developed.  The  kernels 
of  the  integrals  are  bounded  and  well  behaved.  This  type  of 
integral  equation  is  well  understood  in  the  literature  [37,38], 
and  various  means  of  finding  exact  or  approximate  solutions  are 
available. 

Such  a  method  is  developed  in  Chapter  V  for  the  Fredholm 
integral  equations  obtained  from  the  K-L  representation.  From  the 
theory  of  linear  operators,  the  method  of  successive  approximations 
is  applied  to  obtain  an  iterative  power  series  expansion  for  the 
unknowns.  The  expansion  Is  In  terras  of  the  resistivity  of  the  wedge, 
and  converges  uniformly  for  particular  ranges  of  values  of  this 
parameter.  Bounds  for  the  regions  of  a  convergence  are  given.  The 
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chapter  Is  concluded  with  a  discussion  of  how  similar  expansions 
obtained  from  the  difference  equations  may  not  exist  in  certain  cases. 

A  summary  and  discussion  of  the  results  follow  in  Chapter  VI. 
Appendix  A  contains  a  set  of  conditions  for  the  existence  of  the 
representations  used  by  Maliuzhinets  and  Kontorovich  and  Lebedev,  along 
with  their  corresponding  inverses. 

1 .3  Discussion  of  the  Impedance  and  Resistive  Boundary  Conditions 
Before  proceeding  with  an  analysis  of  the  scattering  by  a 
resistive  wedge,  it  is  appropriate  to  review  the  mathematical 
implications  of  the  impedance  and  resistive  boundary  conditions. 

In  vector  form,  the  impedance  boundary  condition  on  the  surface 
S  of  a  body  is  given  by 

E  -  (n  •  E)n  *  nZn  x  R  ;  on  S  ,  (1.1) 


where  (E,H)  are  the  total  fields  in  the  region  surrounding  the 
body,  assumed  to  be  free  space,  Z  is  the  intrinsic  impedance  of  free 
space,  and  n  is  a  unit  vector  normal  to  S  and  directed  into  the  region 
containing  the  fields  (see  Fig.  1.2).  The  dimensionless  parameter  n 
is  the  surface  impedance  of  the  boundary  normalized  to  free  space. 
Physically  speaking,  if  the  body  consists  of  a  material  with  large 
refractive  index,  and  hence  large  relative  complex  permittivity, 
n  has  the  form  [39] 


(1.2) 
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where  e,ji,a  are  the  permittivity,  permeability,  and  conductivity 
of  the  body,  respectively,  and  eQ,  u0  are  the  corresponding  free 
space  parameters.  An  e"1u)t  time  dependence  has  been  assumed  and 
suppressed  here  and  throughout  this  work  with  w  the  angular  frequency 
of  the  EM  field.  Note  that  as  a  -*■  »,  n  0  and  (1.1)  becomes  the 
boundary  condition  for  perfect  conductivity,  i.e., 

n  x  E  =  0  , 

as  expected.  It  was  in  this  context  that  (1.1)  was  introduced  [8] 
as  an  approximation  to  the  boundary  conditions  at  a  body  with 
large  but  finite  conductivity.  A  discussion  of  the  validity  of  the 
approximation  can  be  found  in  [40].  The  utility  of  (1.1)  is  that 
it  reduces  the  problem  of  determining  the  fields  both  inside  and 
outside  the  body  to  that  of  solving  the  exterior  problem  alone, 
subject  to  a  boundary  condition  which  describes  the  material  properties 
of  the  body.  All  interior  fields  are  identically  zero.  As  discussed 
in  Section  1.1,  several  scattering  problems  satisfying  such  a 
condition  have  been  solved. 

Regardless  of  its  physical  implications,  mathematically 
(1.1)  relates  the  tangential  components  of  E  and  fi  via  a  parameter  n. 

In  terras  of  the  equivalent  electric  and  magnetic  surface  currents 

K  =  n  x  H  (1 .3a) 

e 


x  E  * 


(1.3b) 
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respectively,  It  follows  that 

K  =  -nZn  x  K  .  (1.4) 

me 

An  interesting  duality  transformation  exists  for  the  impedance 
boundary  value  problem.  It  can  be  stated  via  the  following  theorem 
[39]: 

Theorem  1 :  If  the  electromagnetic  field  incident  upon  a  body 
satisfying  the  boundary  condition  (1.1)  is  denoted  by 

(t1  .H1 )  ■  (F.YG)  . 

and  the  scattered  field  by 

(is,Hs)  =  (T(n).Yg(n))  . 

where 

(i, H)  *  (E1  +  Es  ,  H1  +  H')  , 

then  for  an  incident  field 

(if  ,H7*)  =  (-ZG.F)  , 

the  scattered  field  is 

(ES,HS)  =  (-Zg(l/n),  f(l/n) )  . 

This  is  equivalent  to  the  transformation  E  -*■  H,  ZH  -►  -YE,  n  -►  l/n» 
and  can  easily  be  derived  by  taking  the  cross  product  of  (1.1)  with 
n,  yielding 


,  -  - 


\ 
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H  -  (n  •  H)n  =  .  1  Yn  x  E  .  (1.5) 

Equation  (1.5)  is  simply  a  restatement  of  (1.1)  under  the  prescribed 
transformation.  By  allowing  n  -*■  =°,  (or  1/n  0),  (1.5)  reduces  to 

n  x  H  *  0  , 

which  is  the  boundary  condition  for  a  perfect  (nonetheless  fictitious) 
magnetic  conductor. 

Suppose  instead  that  the  material  body  in  Fig.  1.2  is  replaced 
by  a  thin  shell  of  thickness  6,  coincident  with  surface  S.  The  shell 
is  composed  of  the  same  material  as  the  body,  with  parameters  e,y,o. 

A  solution  for  the  scattering  from  such  a  structure  requires  the 
determination  of  the  fields  in  the  three  regions  defined  in  Fig.  1.3. 
If,  however,  one  again  considers  the  situation  when  the  complex 
permittivity  becomes  large  while  simultaneously  allowing  the 
thickness  to  decrease,  such  that  their  product  remains  constant,  then 
the  three  region  problem  can  be  approximated  by  a  two  region  problem 
with  an  appropriate  set  of  boundary  conditions  on  the  surface  S. 

A  mathematical  description  of  the  limiting  process  is  given  in  (25,26]. 
Referring  to  the  geometry  in  Fig.  1.4,  the  approximate  vector 
boundary  conditions  that  now  apply  are 
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4  Equivalent  Geometry  for  Fig.  1.3  with  the  Scattering 
Body  Replaced  by  a  Resistive  Surface. 
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n  x  [E]*  =0  (1 .6a) 

n  x  (n  x  E)  =  -Rn  x  [H]*  (1.6b) 

where  the  notation  [  ]*  denotes  the  discontinuity  in  the  value  of  a 
quantity  across  the  surface  S,  or  more  exactly, 

[f]*  *  -  f2  ;  on  S  , 

with  f  an  arbitrary  function.  The  subscripts  refer  to  the 
corresponding  values  in  Regions  1  and  2  of  Fig.  1.4.  In  terms  of 
equivalent  surface  currents 


Re  *  n  x  [H]*  (1.7a) 

iS,  -  -n  x  [E]*  ,  (1.7b) 

(1.6)  becomes 

K,,,  -  0  (1.8a) 

n  x  (n  x  E)  *  -RK  .  (1 .8b) 

e 


The  parameter  R,  the  resistivity  of  the  layer,  is  given  by 
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in  the  limit  as  described  above.  The  resistivity  R  has  units  of 
ohms  per  square  in  the  MKS  system.  When  a  »  e,  R  tends  toward  a 
purely  real  number,  given  by  [263 


«  -  h  ■  n.io) 

In  general,  for  a  passive  material  with  a  >_  0,  e  >_  eQ,  R  is  a 
complex  number  lying  in  the  first  quadrant  of  the  complex  plane.  From 
(1.7),  a  resistive  surface  is  equivalent  to  an  electric  current  layer 
whose  strength  is  proportional  to  the  tangential  electric  field 
at  its  surface.  Since  the  tangential  electric  field  is  continuous 
across  the  layer  {from  (1.6a)),  the**e  are  no  magnetic  currents. 

It  is  interesting  to  note  that  in  the  limit  as  R+  0,  (1.6)  becomes 
a  perfectly  conducting  boundary  condition,  while  for  R  ■+>  ®,  (1.6)  can 
be  written  as 

n  x  [E]*  *  0 
n  x  [H]*  =  0  . 

These  continuity  conditions  are  satisfied  by  the  incident  field 
alone,  and  hence  there  is  no  scattered  field.  Equivalently,  the 
scattering  body  has  ceased  to  exist. 

As  was  done  for  the  impedance  boundary  condition  (1.1),  It  Is 
convenient  to  normalize  the  resistivity  via  a  dimensionless  parameter 


0.11) 


n 


2R 
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The  factor  of  two  Is  introduced  because  of  an  interesting  result 
from  half-plane  diffraction  [27].  The  total  electric  current  on  a 
half  plane  with  surface  impedance  n  is  identical  to  that  on  a  resistive 
half-plane  with  R  satisfying  (1.11).  However,  this  result  cannot  be 
generalized  for  arbitrary  geometries. 

Recalling  the  duality  transformation  presented  in  Theorem  1, 
its  application  to  (1.6)  yields 

n  x  [H3t  =  0  (1.12a) 

n  x  (n  x  H)  =  R*n  x  [E]*  ,  (1.12b) 

where 

R*  »  3j-  •  (1.13) 

While  (1.5)  is  simply  a  restatement  of  the  impedance  boundary  condition 
under  the  prescribed  duality  transform,  i.e.,  an  impedance  boundary 
condition  is  its  own  dual,  examination  of  (1.12)  indicates  a  similar 
analogy  does  not  exist  for  the  resistive  boundary  condition.  Indeed, 
(1.12)  describes  a  surface  characterized  by  a  jump  discontinuity  in 
the  tangential  electric  field,  but  no  discontinuity  in  the  tangential 
magnetic  field.  Such  a  surface  supports  only  a  magnetic  surface 
current  1^,  and  has  been  referred  to  as  a  "conductive"  sheet,  with 
conductivity  R*  mhos  per  square  [27].  A  discussion  of  the  duality 
relations  for  Impedance  and  resistive  boundary  conditions  Is 
presented  by  Senior  [41]  in  relation  to  Babinet's  principle. 
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It  Is  appropriate  to  further  compare  the  impedance  and 
resistive  boundary  conditions  In  the  context  the  problem  considered 
here.  A  wedge  with  upper  and  lower  faces  S+  and  S_  respectively  is 
illuminated  by  a  plane  electromagnetic  wave  incident  normal  to  the 
edge  of  the  wedge.  By  considering  the  cases  in  which  the  incident 
electric  or  magnetic  field  is  parallel  to  the  edge  (E-  or  H-polarization 
respectively),  as  illustrated  in  Figs.  1.5a  and  1.5b,  the  vector 
scattering  problem  can  be  reduced  to  two  scalar,  two-dimensional 
scattering  problems. 

In  the  event  that  S+  and  S_  satisfy  impedance  boundary  conditions, 
(1.5)  for  E-polarization  becomes 

l^  +  ^-U  =  0  ;  on  S+,S_  (1.14a) 

where  E  =  uU,  u  being  a  unit  vector  parallel  to  the  edge  of  the 
wedge.  The  normal  derivative  3/an  is  defined  as 


i.e.,  it  is  taken  in  the  direction  of  the  unit  vector  n,  which  from 
Fig.  1.4  is  directed  toward  Region  1  containing  the  incident  field. 
The  propagation  constant  k  is  defined  as 

k  8  c  ’ 

where  c  is  the  velocity  of  light  in  vacuo. 
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Obviously,  (1.15)  and  (1.16)  are  not  duals  of  each  other,  and 
therefore  the  scattering  problem  for  an  arbitrary  polarized  incident 
plane  wave  requires  the  solution  of  two  different  scalar  problems. 
Further  complications  arise  from  the  existence  of  fields  in  both 
Regions  1  and  2,  as  implied  by  the  bracketed  quantities  in  (1.15) 
and  (1.16). 

For  a  more  detailed  comparison  of  the  impedance  and  resistive, 
as  well  as  "conductive1'  boundary  conditions,  the  reader  is  referred 
to  the  works  of  Senior  [27,29,30,41,42], 


CHAPTER  II.  FORMULATION  OF  THE  SCATTERING  BY  A  RESISTIVE  WEDGE 


2.1  Statement  of  the  Problem 

Having  provided  a  qualitative  analysis  of  the  resistive 
boundary  condition,  including  its  physical  and  mathematical  implications, 
and  having  compared  it  to  more  conventional  boundary  conditions,  a 
rigorous  formulation  of  the  electromagnetic  scattering  of  a  plane  wave 
normally  incident  on  a  resistive  wedge  is  presented. 

The  geometry  under  consideration  is  shown  in  Fig.  2.1.  A  wedge 
of  included  angle  2\p  composed  of  two  resistive  sheets  S+  and  S_  has  its 
vertex  coincident  with  the  z-axis  of  a  cylindrical  coordinate  system 
(p,<t>,z).  The  resistivity  R  of  the  sheets  is  a  scalar  constant.  The 
azimuthal  coordinate  <j>  is  chosen  to  take  on  the  values  -iji  <  $  <  2ir  -  \p . 
The  upper  face  of  the  wedge  lies  in  the  plane  <f>  =  ^;  the  lower  face  in 
the  plane  <p  =  2n  -  ip.  The  exterior  region  (Region  1)  is  defined  by 


(p,4>,z)  C  Region  1 


fO  <  P  <  » 

ip  <_  <_  2ir  -  ip 

-00  <  Z  <  00 


Similarly,  the  interior  region  (Region  2)  is  defined  by 


(p,$,z)e  Region  2  =^> 


[0  <  P  <  - 

“>/'<_  <p  <_  ip 

<  z  <  « 
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Fig.  2.1  Resistive  Wedge  of  Included  Angle  2\p  Illuminated  by 
a  Plane  Wave  at  Normal  Incidence. 
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Region  1  and  2  are  assumed  to  be  free  space,  \**ith  e  =  eQ, 
u  =  pq.  An  E(H)-  polarized  electromagnetic  plane  wave  with 


is  normally  incident  upon  the  wedge  from  an  angle  $  ,  which  lies  in 
Region  1  (see  Fig.  2.1).  Because  of  the  symmetry  of  the  geometry 
about  the  plane  y  =  0,  it  is  sufficient  to  consider,  without  loss  of 
generality,  the  following  relation  between  the  angle  of  incidence 
and  the  half  angle  of  the  wedge: 

0  <  Ip  <  <j>  <  TT 

0  — 

Since  the  entire  problem  is  independent  of  the  z  coordinate,  the 
two-dimensional  problem  in  the  plane  z  =  0  will  be  considered  from 
this  point  on.  As  stated  previously,  a  harmonic  time  variation  of 
the  form  e’1ut  is  assumed  and  suppressed  throughout. 

Because  the  problem  is  two-dimensional,  the  total  electro¬ 
magnetic  field  (E,H)  can  be  determined  from  the  two  scalar  quantities 
Ez  and  Hz-  By  assuming  in  turn  the  E-  and  H-polarized  incident 
fields  (2.1),  (for  which  Hz  and  Ez  are  zero,  respectively),  the 
solution  to  the  general  problem  is  reduced  to  solving  two  scalar 
problems  for  a  single  unknown. 

As  implied  in  (2.1),  the  unknowns  Ez,Hz  will  be  denoted 


by 
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Ez  =  U  =  U1  +  Us  (2.2a) 

H2  =  V  =  V1  +  Vs  (2.2b) 

where  the  superscripts  i  and  s  denote  the  incident  and  scattered 
fields,  respectively.  The  fields  U  and  y  are  functions  of  the 
variables  (p,<j>),  and  as  a  result  it  is  convenient  to  distinguish 
their  values  in  Regions  1  and  2  via  appropriate  subscripts,  viz, 

U  5  Uj  s  uj  +  ,  for  (p,ij>)£  Region  1  (2.3a) 

=  U2  =  ul  +  uf  ,  for  (p,<t>)E  Region  2  ,  (2.3b) 

and  similarly  for  V.  As  is  traditionally  the  case,  the  incident  field 
is  assumed  to  permeate  the  entire  space. 

From  Maxwell's  curl  equations  in  free  space 


7  x  E  3  iuvijjH  , 
V  x  H  =  - iwEgl  • 


the  following  results  can  be  shown  in  two  dimensions  (3/3z  =  0): 
E-polarization:  Ep  =  =  Hz  =  0  ,  Ez  =  U 


ikZH 

p 

ikZH 

9 


=  1  5U 

p  H 

3p 


(2.4a) 

(2.4b) 
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H-polarization:  H  =  H.  =  E_  =  0  ,  H_  =  V 

— c -  p  #  z  z 

(2.5a) 
(2.5b) 


=  _1  3V 

P  3$ 
3V 

3p  ' 


In  addition,  it  can  be  shown  that  both  U  and  V  satisfy  the  scalar 
Helmholtz  equation 


(2.4) 


in  Regions  1  and  2,  where  v2  is  the  two-dimensional  scalar  Laplacian. 
Expressing  the  resistive  boundary  conditions  (1.15)  and  (1.16)  in 
polar  coordinates  (p,<j>),  the  following  conditions  on  the  surfaces  S+ 
and  S_  hold  for  E-polarization: 


Similarly,  for  H-polarization  the  conditions  are 
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(2.7) 


Because  the  medium  and  the  scattering  body  are  both  infinite 
in  extent,  the  functions  U  and  V  are  required  to  satisfy  an 
additional  constraint  as  p  »  in  the  form  of  a  Sommerfeld  radiation 
condition.  Specifically,  by  decomposing  U  and  V  into  a  geometrical 
optics  component,  obtained  from  ordinary  ray  theory,  and  a  diffracted 
component,  viz 

U  =  U9  +  Ud  (2.8a) 

V  =  V9  +  Vd  (2.8b) 


each  of  which  is  a  discontinuous  function  of  <p ,  then  the  diffracted 
components  satisfy* 


★ 


The  reasoning  behind  this  decomposition  is  discussed  by  Williams  [11]. 
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uniforraly  for  all  4>  appropriate  to  Regions  1  and  2.  Condition  (2.9) 
is  discussed  in  [6]  and  elsewhere;  its  vector  analog,  the  Silver- 
Muller  condition,  can  be  found  in  several  texts  on  electromagnetic 
theory,  see  Jones  [43],  for  example.  Roughly  speaking  (2.9)  is 
equivalent  to  requiring  that  the  diffracted  field  have  the  form 
of  an  outgoing  cylindrical  wave  at  infinity,  decaying  as  the 
reciprocal  of  the  square  root  of  the  distance  from  the  line  z  =  0. 

Finally,  the  geometrical  singularity  presented  by  the  vertex 
of  the  wedge  requires  the  specification  of  an  additional  condition 
governing  the  behavior  of  the  fields  in  the  vicinity  of  the  edge. 
This  physical  constraint,  known  as  the  edge  condition,  is  usually 
expressed  by  requiring  that  the  stored  electric  and  magnetic  energy 
in  any  neighborhood  of  the  edge  be  finite;  that  is 


u0i£[2  +  p0l"i2>dv  =  0 


(2.10) 


Jones  [43]  derives  a  uniqueness  theorem  for  finite  dielectric  bodies 
and  infinite  perfectly  conducting  bodies  based  on  (2.9)  and  (2.10).  He 
indicates,  however,  that  these  equations  do  not  appear  to  be 
sufficient  to  insure  uniqueness  for  infinite  dielectric  structures 
(such  as  a  dielectric  wedge)  [44].  While  the  uniqueness  of  the 


resistive  wedge  problem  has  not  been  addressed  in  the  literature, 
it  will  he  assumed  that  a  solution  satisfying  (2.1)  through  (2.10) 
provides  the  desired  results.  This  assumption  is  legitimately 
questioned,  since,  as  was  pointed  out  in  Chapter  I,  there  exist 
similarities  between  the  resistive  and  dielectric  wedge  scattering 
problems. 

There  are  several  physically  intuitive  means  of  expressing 
the  edge  condition  [45,46]  which  are  equivalent  to  (2.10)  for  a 
perfect  conductor.  These  include  zero  induced  charge  or  finite 
surface  current  and  charge  conditions.  An  overview  of  these  works 
and  a  general  discussion  of  the  edge  condition  for  perfect 
conductors  can  be  found  in  Oones  [43], 

It  is  relatively  easy  to  show  for  a  straight  edge  that  (2.10) 
is  equivalent  to  requiring  that  no  component  of  the  fields  be  more 
singular  than  p~5  in  the  neighborhood  of  the  edge,  where  6  <  1. 

More  exactly, 


max  { f.|}  =  0(p"5)  ;  p  -*■  0  ,  (2.11) 

where  { f ^ >  are  the  components  of  the  electric  and  magnetic  fields,  and 
0(  )  is  the  standard  order  relation,  whereby 


p  -*■  0 

lim  ( fg"1 ]  <  A  . 
p  -►  o  '  ' 


f  =  0(g)  ; 

3  A  *  constant 
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Meixner  and  others  [47,48]  have  investigated  the  possible  values 
of  the  parameter  d  for  various  wedge-like  regions,  including  both 
perfectly  conducting  and  dielectric  wedges.  By  assuming  series 
expansions  in  powers  of  p  for  the  various  field  components,  Meixner 
derives  expressions  for  the  series  coefficients  and  admissible  values 
of  5  via  Maxwell's  equations  and  the  boundary  conditions.  In 
particular,  he  shows  that  the  lowest  order  (most  singular)  terms  are 
equivalent  to  those  obtained  for  the  corresponding  static  field 
problem.  An  important  consequence  of  this  result  is  that  the  field 
components  parallel  to  the  edge  (in  this  case,  Ez  and  Hz)  are  finite 
for  all  geometries.  Although  some  doubt  has  recently  been  expressed 
regarding  the  validity  of  Meixner's  series  expansion  in  the  dynamic 
(io  f  0)  case  [49],  it  appears  that  the  lowest  order  terms  are  still 
the  correct  asymptotic  forms. 

Thus,  for  the  resistive  wedge,  it  is  assumed  that 

U,V  =  0(p6)  ;  6  >  0  ,  p  -  0  ,  (2.12) 

uniformly  for  all  $  in  each  of  Regions  1  and  2.  Note  that  in  the 
event  that  6  =0,  i.e.,  U,V  =  constant  at  p  =  0,  equation  (2.6) 
guarantees  that  the  value  of  the  constant  is  the  same  in  Regions  1 
and  2  for  E-polarization.  However,  the  same  is  not  true  for 
H-polarization  (see  (2.7));  there  may  be  a  discontinuity  in  the 
value  of  V  at  p  =  0  in  passing  from  Region  1  to  Region  2.  This 
observation  is  related  to  the  fact  that  a  resistive  layer  can 
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support  only  an  electric  current,  which  is  radially  directed  for 
H-polarization.  In  the  event  that  this  current  is  zero  at  the 
edge,  then  the  discontinuity  disappears  and  the  value  of  V  at  p  *  0 
is  the  same  in  both  regions. 

In  summary,  the  problem  addressed  in  this  work  may  be  stated 
as  follows.  For  the  geometry  of  Fig.  2.1,  illuminated  by  alternately 
E-  and  H-polarized  plane  waves  of  the  form  (2.1),  two  scalar 
unknowns  U  and  V  are  sought  which  satisfy  the  Helmholtz  equation 
(2.4),  the  boundary  conditions  (2.6)  and  (2.7),  respectively,  and 
the  radiation  and  edge  conditions  (2.9)  and  (2.10),  or  equivalently, 
(2.12).  The  problem  will  be  analyzed  by  applying  the  function- 
theoretic  techniques  of  Maliuzhinets  and  Kontorovich  and  Lebedev. 

It  is  worthwhile  to  note  that  the  more  conventional  method  of  mode 
matching,  otherwise  known  as  the  method  of  separation  of  variables, 
is  not  applicable  to  resistive  (or  impedance)  wedge  scattering 
problems  unless  a  resistivity  (or  impedance)  which  varies  in  a  specified 
manner  as  a  function  of  p  is  assumed  (see  [50],  for  example).  Since 
a  constant  resistivity  is  assumed  in  this  work,  the  mode  matching 
method  is  not  a  valid  approach  to  the  problem. 

Before  proceeding  with  the  analysis,  a  decomposition  for  the 
unknowns  will  be  outlined  which  will  prove  instrumental  in  simplifying 
several  of  the  resulting  expressions,  especially  the  boundary  conditions 
(2.6)  and  (2.7). 
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2.2  The  Method  of  Symmetric  and  Antisymmetric  Components 

It  is  a  well  known  result  that  an  arbitrary  function  f(x) 
defined  on  an  interval  (a,b)  can  be  decomposed  into  two  functions 
which  are  respectively  even  (symmetric)  and  odd  (antisymmetric)  about 
the  midpoint  c  of  (a,b).  Denoting  these  functions  by  fe  and  f°,  and 
letting 
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An  important  consequence  of  this  result  is  that  a  knowledge  of  the 
single  function  f(x)  on  (a,b)  is  equivalent  to  a  knowledge  of  two 
functions,  fe(x)  and  f°(x),  on  the  half  interval  (a,cj. 

The  exterior  and  interior  regions  of  the  resistive  wedge 
(Fig.  2.1)  can  be  considered  as  intervals  [(f,2TT-if]  and  in  the 

angular  variable  <f.  Likewise,  the  unknowns  and  ll2,V2  in 

these  regions  represent  functions  of  $  on  the  appropriate  intervals. 
Hence,  the  decomposition  (2.13)  can  be  applied  to  the  unknowns  in 
Regions  1  and  2  by  writing 

U  (p,<f)  =  i  [Ue  (p,*)+U°  (p ,<*>)]  (2.17a) 

1  »2  ‘  1  »2  1  »2 

V  (Pi*)  =  i  [Ve  (p  »<f> )  +  V°  (p,<f)]  ,  (2.17b) 

1  *2  ^  1.2  1.2 

where,  neglecting  the  p-dependence, 

U*’°U)  =  U  (♦)  +  U  (2ir  -  <f)  (2.18a) 

Ue’°(<f)  =  U  (<f)  !  U  (-(f)  (2.18b) 

2  2  2 

and  similarly  for  V  and  V  .  The  plus  sign  corresponds  to  the 
1  2 

symmetric  (e)  component,  and  the  minus  sign  to  the  antisymmetric  (o) 

component.  From  (2.15)  it  is  obvious  that  Ue,Ve  (U°,V°)  are  even  (odd) 

1111 

about  <f  =  iti  while  Ue,Ve  (U°,V°)  are  even  (odd)  about  f  =  0.  Of 

2  2  2  2 

greater  significance  is  equation  (2.16),  from  which  It  is  determined 
that 
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and 


It  (ue.ve) 

3<J>  11 


=  0 


4>  — 7T 


(U°,V°) 
1  1 


*  0 


4>si t 


~  (ue’v6) 

3<j>  2  2 


=  0 


4>  —o 


(U°,V°) 

2  2 


$*0 


(2.19a) 

(2.19b) 

(2.20a) 

(2.20b) 


It  is  beneficial  to  recall  equation  (2.3),  in  which  U  and  V 
were  expressed  as  the  sum  of  an  incident  and  scattered  field.  The 
same  result  holds  for  Ue,  U°,  Ve,  and  V°,  i.e.. 


-  - 

- 

ue 

uei 

+  ues 

1  >2 

1  >2 

1  »2 

u° 

U°i 

+  uos 

1  »2 

1  >2 

1  >2 

(2.21) 


and  likewise  for  V®  2>  V°  2-  It  is  hoped  that  the  reader  will  forgive 

the  rather  cumbersome  use  of  superscripts  and  subscripts  which 

the  author  has  employed.  From  (2.18)  and  (2.1)  it  is  simple  to  show 


that 
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yei  =  yei  =  yei  s  yei  =  eikpcosU-4>0)  +  ei  kpcos  (4»-N>0)  C2.22a) 

uOT  =  yOi  a  yoi  3  yoi  =  eikpcos(*-*0)  _  eikPcos(<f>+4>0)  (2.22b) 

Therefore  the  notation  for  the  various  incident  field  components 
can  be  simplified  by  denoting  the  right-hand  sides  of  equations 
(2.22a,b)  by  Ie  and  1°,  respectively.  Equation  (2.21)  now  has 
the  less  confusing  form 


r 

—  — 

Ul,2 

Ie  +  ues 

1  ,2 

u?  , 

1°  +  u?s„ 

1,2 

- 

1,2 

ve 

Vv 

Ie  +  vf2 
*  * 

- 1 

CM 

•* 

o  ^ 

> 

- 1 

T°  +  V0S  * 

l.2  V 

(2.23) 


From  the  symmetry  of  the  geometry  it  is  easy  to  show  that  the 
left-hand  side  of  (2.23)  satisfies  the  boundary  conditions  of  the 
original  problem  (2.6)  and  (2.7),  as  well  as  the  wave  equation  and 
radiation  and  edge  conditions.  Thus  (U,  ,,V.  „)  are  the  E-  and  H- 
polarized  solutions  to  the  scattering  problem  with  incident  field  Ie, 
while  (U°  2,V°  2)  are  the  E-  and  H-polarized  solutions  for  the 
incident  field  1°.  Furthermore,  by  essentially  replacing  the  boundary 
conditions  on  the  lower  wedge  surface  S_  by  the  conditions  (2.19)  and 
(2.20),  it  is  sufficient  to  determine  the  unknowns  in  the  upper  half 
space  y  >_  0  to  completely  solve  the  problem,  as  illustrated  in 
Fig.  2.2. 
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*+  / 

u’t*)  ' 

U’U) 

^  / 

<  =* 
^  \ 

\ 

u1  (2n-*)  s  \ 

/ 

/ 

/ 


/ 


-A 


/ 


m  *  0 

3<j> 


(a) 


/ 


U1  (4>) 


-U1  (2tt-4>  > 


/ 


< 


\ 


\ 


u1^)  / 

^  / 

r  >  >  )  i  >  rr  t  >  i  >  >  rrrf 


U  =  0 


(b) 

».•  »=  asss  i^nniseT^s?aibt'«^c^api«  ■&?-.. 

Analogous  Results  hold  for  V($). 


! 

i 


r 


I 
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The  new 


boundary  conditions  for  U®>2  become 


3U®  3U®' 

1  + 

34~  “ 

L  J 

4=4 

4*4 


w;  -  uj] 


«  0 


4=4 


5 

3<t > 


3U' 

If 


» 


4>*o 


=  0  , 


(2.24a) 

(2.24b) 

(2.24c) 

(2.24d) 


while  for  V®  ,  the  following  hold: 
1  »2 


1  as 

P  34 


1 

I 

I 

I 

I 

I 


4=4 

1_ 
p 


+  ^  [V*  ■ 

■  V®]  ■  0 

z  4*4 

(2.25a) 

1  “  ° 

(2.25b) 

34  34 

‘4*4 

—1 

*  0 

(2.25c) 

3*  w 

uSi 

a  0  . 

(2.25d) 

**  U-o 

-38- 


For  U°  2,V°  2,  equations  (2.24c,d)  and  (2.25c,d)  are  respectively 
replaced  by  the  primed  equations 


U 


o 

1 


4>=TT 


=  0 


(2.24c') 


and 


(2.24d') 

(2. 25c') 


V° 


(j)=0 


=  0 


( 2 . 25d  * ) 


In  the  remaining  equations  the  symmetric  (e)  components  are  simply 
replaced  by  the  antisymmetric  (o)  components. 

The  most  important  consequence  of  the  modified  conditions 
(2.24  through  2.25)  is  the  presence  of  only  a  single  "mixed" 
transition  equation  for  each  unknown,  i.e.,  equations  (2.24a)  and 
(2.25a).  These  conditions,  which  contain  the  only  dependence  on 
the  resistivity  n»  are  most  responsible  for  the  complexity  of  the 
problem.  The  method  of  symmetric  and  antisymmetric  components  has 
essentially  eliminated  this  equation  for  the  lower  sheet.  The 
implications  of  this  result  will  become  apparent  in  Chapters  III 
and  IV,  in  which  the  scattering  problem  will  be  analyzed  via 
Maliuzhinets'  method  and  the  Kontorovich-Lebedev  transform. 


CHAPTER  III.  THE  METHOD  OF  MALIUZHINETS 


By  far  the  most  successful  technique  for  solving  boundary 
value  problems  in  wedge- like  regions  was  put  forth  by  G.  D. 
Maliuzhinets  in  his  doctoral  thesis  [13]  and  subsequent  publications 
[14,15].  His  method,  which  fundamentally  is  a  generalization  of  the 
method  of  images  described  by  Carslaw  [3],  was  able  to  provide  the 
solution  for  the  scattering  from  an  impedance  wedge  with  differing 
impedances  on  each  of  the  faces  [15].  Various  aspects  of 
Maliuzhinets1  technique  were  inadvertently  and  independently  described 
by  Senior  [16]  and  Williams  [17]  in  solving  similar  problems. 

The  basis  for  Maliuzhinets'  method  is  the  ability  to  represent 
the  total  scalar  fields  (U,V)  as  a  Sommerfeld  integral  of  the  form 


where  y  Is  a  contour  in  the  complex  a-plane,  consisting  of  two  loops 
Y  and  y2»  symmetric  about  the  point  a  =  it,  as  Illustrated  in  Fig.  3.1. 
A  contour  integral  of  this  type  was  first  employed  by  Sommerfeld  in 
his  classic  solution  for  the  scattering  of  a  plane  wave  by  a  perfectly 
conducting  half-plane  [1],  and  as  a  result  (3.1)  is  often  referred 
to  as  a  Sommerfeld  Integral.  Maliuzhinets  has  shown  that  (3.1)  has 
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a  unique  inversion  (s,t)  within  a  particular  set  of  functions, 
provided  (U,V)  satisfy  certain  boundedness  conditions  at  p  »  0  and 
p  ->■  ».  An  outline  of  Maliuzhinets1  uniqueness  theorem  is  given  in 
Appendix  A. 

A  direct  consequence  of  this  uniqueness  theorem  is  that  a 
nonzero  solution  of  the  homogeneous  equation 

f(a)eikpC0Sa  da  =  0  (3.2) 

must  be  even  about  the  point  a  =  it,  i.e.. 


f(a)  =  f (2tt  -  a) 


(3.3) 


provided  f(a)  =  0(e^"a^  Inta^ )  as  ]Im  a|  -►  »  within  the  loops  Yj  and 
y  ,  where  a  >  0.  This  result  is  also  derived  in  Appendix  A.  The 
order  relation  above  holds  for  functions  of  (p,<j>)  satisfying  the 
radiation  and  edge  conditions  (2.9)  and  (2.10),  respectively. 

Recalling  the  representation  (3.1)  for  (U.V),  it  is  easy  to 
show  that  indeed  (U.V)  satisfy  the  Helmholtz  equation 


(v2  +  k2) 


(3.4) 


In  order  to  apply  (3.1)  to  the  boundary  conditions  (2.6)  and  (2.7) 
(note  the  (3.1)  is  a  representation  for  the  total  fields  ( U , V ) ) , 
it  is  appropriate  to  first  consider  the  normal  derivative 


s(a  -  i p) 

JkpCOSa  3_  da 

**  t(a  —  <j> ) 


ikpcosa  3 


J"  sta  -  <(>) 


t(a  -  <|>) 


(3.4) 


The  interchange  of  the  order  of  differentiation  and  integration  is 
justified  by  the  uniform  convergence  of  the  integral  along  y  provided 
(s,t)  satisfy  the  order  relation  mentioned  previously.  Integration 


by  parts  leads  to 


1  i_ 

P  3<J> 


1 _  Oikpcosa 

IPfF  e 


s(a  -  $) 


t(a  -  $) 


«  -*■  ya 


sin  a  e 


ikpcosa 


S (a  -  <p) 


t(a  -  4> ) 


(3.5) 


where  denotes  the  various  endpoints  of  the  contour  y  as 
] Im  a|  •+  «.  From  the  behavior  of  (s,t)  and  the  exponential  (for 
p  >  0)  within  the  shaded  portions  of  Fig.  3.1,  it  can  be  shown  that 
the  first  terra  of  equation  (3.5)  is  zero,  and  thus 


1  £_ 

P  <*<t> 


.  sin  a  eikpC0Sa 


s(a  -  <J>) 


da  .  (3.6) 


t(a  - 
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Substitution  of  (3.1)  and  (3.6)  into  the  boundary  conditions  (2.6)  and 
(2.7)  results  in  the  following  equations  for  (s,t):  for  E-polarization 


+  -  s  (a  -  Z*  + 
n  i 


J  £  S^(a  -  2ir  +  ip)  -  S^ (a  +  ^)  j  e 


ikpcosa 


da 


and  for  H-polarization: 


J  ^  sin  at  (a  -  *)  -  |-|ti(a  -  ♦)  -  t2(o  -  ♦)  j  j 

J  s 1 n  a  £  t^ (a  -  if»)  -  t^(a  -  ^  )  j  e 

Y 

j  £sin  a  Ma  -  2tt  +  ^ )  +  J  [^(a  -  2»  +  i|>) 


Ikpcosa 


eikpCOSa  do  .  0 

*  0 

(3.7) 


eikpCOSa  ^  =  0 

da  »  0 


-  t  (a 

2 

t(a-2rr+i|>)-t(a  + 
1  2 


♦  *))  ]  excuse 

♦  )]  »1kt’cosa  * 


*  0 


0 

(3.8) 


In  equations  (3.7)  and  (3.8)  the  subscripts  1,2  on  (s.t) 


correspond  to  the  values  of  the  functions  in  Regions  1,2, 
respectively. 

Following  Maliuzhinets'  derivation  for  the  impedance  wedge  [15], 
and  by  virtue  of  equations  (3.2),  (3.3),  and  the  results  of  Appendix 
A,  equations  (3.7)  and  (3.8)  are  equivalent  to  the  following  set  of 
functional  equations: 

E-polarization: 

n  Sin  a[s  (a  -  i|i)  -  S  (a  -  ip ) 3  -  2s  (a  -  ij>) 

1  2  1 

=  -n  sin  a[s  (2ir  -  a  -  <|/)  -  S  (2ir  -  a  -  ip ) ]  -  2s  (2ir  -  a  -  ip)  (3.9a) 

1  2  1 

S  (ct  -  4>)  “  S  (a  -  )  *  S  (2it  -  a  -  <p )  -  S  (2tt  -  a  -  tp ) ( 3 . 9b) 

12  1  2 

n  sin  a[s  (a  -  2n  +  ij>)  -  s  (a  +  ip ) 3  +  2s  (a  -  2ir  +  i|>) 

1  2  1 

*  -n  sin  a[s  (-a  +  ty)  -  s  (2w  -  a  +  i(»)]  +  2s  (-a  +  ip)  (3.9c) 
12  1 

S  (a  -  2ir  +  ip )  -  S  (a  +  ip )  *  S  (-a  +  tp )  -  S  (2*  -  a  +  ip ) ( 3 . 9d ) 

1  2  2  2 

H-polarization: 

2  sin  a  t  (a  -  xp)  -  n[t  (a  -  i|»)  -  t  (a  -  \|»)] 

1  1  2 

*  -2  sin  a  t  ( 2tt  —  a  —  ip )  —  n[t  (2tt  -  a  *  ip)  -  t  ( 2ir  -  a  -  ip ) 3  (3.10a) 
1  1  2 
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t  (a  -  ip)  -  t  (a  -  \l>)  -  -t  (2tt  -  a  -  0)  +  t  (2rr  -  a  -  i/O  (3.10b) 

12  1  2 

2  sin  a  t  (a  -  2tt  +  ij/)  +  n[t  (a  -  2ir  +  </0  -  t  (a  +  i(0J 
1  1  2 

*  -2  sin  ct  t  (-a  +  ij>)  +  n[t  (-a  +  i(0  -  t  (2ir  -  a  +  iji)]  (3.10c) 
1  1  2 

t  (a  -  2ir  +  ip)  -  t  (a  +  <|i)  =  -t  (-a  +  <(0  +  t  (2ir  -  a  +  t|»)  (3.1  Od ) 

1  2  12 

Equations  (3.9)  and  (3.10)  can  be  put  in  simpler  terms.  Adding 
(3.9a)  and  (3.9b)  and  subtracting  (3.9c)  and  (3.9d)  results  In  two 
functional  equations  of  similar  form: 

(1  -  n  sin  a)S  (a)  =  (1  +  n  sin  a)S  (2tt  -  a)  (3.11a) 

1  l 

(1  -  n  sin  o)S  (a)  =  (1  +  n  sin  a)s  (2tt  -  a)  ,  (3.11b) 

2  2 

where 


S  (a)  =  s  (a  -  ty)  -  S  (2tt  -  a  -  i/O  (3.12) 

1  1  2 

S  (a)  *  s  (-a  +  i(0  -  S  (a  +  i(0  .  (3.12) 

2  1  2 

Analogously,  equations  (3.10a)  through  (3.1 Od )  can  be  reduced  to 


(n  -  sin  a)T  (a)  *  (n  +  sin  a)T  (2ir  -  a) 

l  l 


(3.13a) 


(n  -  sin  a)T  (a)  *  (n  +  Sin  a)T  (2ir  -  a) 

2  2 


(3.13b) 


«-»■ 
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where 


T  (a) 

=  t  (a  -  tjj)  +  t  ( 2tt  -  a  -  ij») 

1 

1  2 

T  (a) 

=  t(-a+^)+t(a+ij>) 

(3.14) 

2  1  2 


Hence,  equations  (3.11)  and  (3.13)  represent  a  set  of  coupled  first- 

order  functional  equations  for  the  unknowns  s  ,s  and  t  ,t 

12  12 

respectively.  An  important  observation  is  that  $  »S  ,T  ,T 
satisfy  identical  equations,  and  therefore  can  differ  at  most  by  a 
multiplicative  factor  M(a)  satisfying  (see  Chapter  2  of  [36]) 

M(a)  =  M(2ir  -  a)  .  (3.15) 


Having  determined  S  and  S  ,  one  can  obtain  an  uncoupled 

1  2 

set  of  first-order,  inhomogeneous  functional  equations  for  and 

s  ,  viz 

2 


S^(a  -  i/i)  -  (a  -  2v  +  3ip )  * 

S  (a  -  iji)  -  S  (a  -  2it  +  3^)  = 

2  2 

A  similar  result  holds  for  t  and  t  , 

1  2 

t  (a  -  \p)  -  t^  (o  -  2ir  +  3iJj)  = 

t  (a  “  i|i)  -  t  (a  -  2ir  +  3i^)  = 


S  (a)  -  S  (2ir  -  a  -  2ip )  (3.16a) 

1  2 

S  (a  -  2ir  +  2iji)  -  S  ( 2ir  -  a) 

2  1 

(3.16b) 

i.e., 

T(a)-T(2ir-a-2ifi)  (3.17a) 
1  2 

-T  (a  -  2rr  +  2tfi)  +  T  ( 2ir  -  a) 

2  1 


(3.17b) 
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Once  again  there  is  a  similarity  among  the  functional  equations 

(3.16)  and  (3.17)  describing  s  and  t  .  Indeed,  if  x(a)  is 

1,2  1,2 

an  arbitrary  solution  of  the  homogeneous  equation 

x(ct  -  *)  -  x(<*  -  2tr  +  3*)  =  0  (3.18) 

then  s1>S2»t1>t2  are  of  t^5e  form 


x(a)  +  P(a)  , 


(3.19) 


where  P (a)  is  a  solution  of  the  appropriate  inhomogeneous  equation 

(3.16)  or  (3.17).  That  is  to  say,  s  and  t  are  determined 

1 >2  1 >2 

up  to  an  additive  function  of  period  2n  -  4\p . 

In  order  to  uniquely  determine  s,  „,t  „  additional  constraints 

1*2  1  »2 

are  required.  One  such  constraint  is  the  order  relation  previously 
Imposed: 


1.2 


1,2 


0(e(l-a)|Im  a|)  . 


Im  a  • 


(3.20) 


within  and  on  the  loops  y  and  y  .  Furthermore,  as  is  shown  in 

1  2 

Appendix  A,  for  an  incident  field  of  the  form 


U1  =  y1  =  e~  1  kpcos (<|>  -  *Q)  ^  {3.2}) 
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the  functions 


s.  Cot  -  it) 

-  (ot  -  2tt  +  $  )“ 

tj(a  -  it)  j 

must  be  regular  in  the  strip  <_  Re  a  <_  2tt  -  while 


(3.22a) 


S2(a) 

t2(a) 


(3.22b) 


must  be  regular  in  the  strip  ir  -  ^  Re  at  <_  tt  +  ^ .  Conditions  (3.20)  and 
(3.22)  serve  to  limit  the  choices  of  xM  and  P(o),  as  well  as  M(a) 

(see  (3.15))  in  determining  the  unknowns  SjjS-.tj.t,. 

In  the  event  that  it  may  be  advantageous  to  preclude  the 
derivation  of  the  corresponding  conditions  on  Sj  .Sj.Tj.Tj,  as  well  as 
to  avoid  the  large  number  of  arbitrary  unknowns  x,P,M  and  so  on,  it 
is  possible  to  eliminate  S2  and  t2  from  equations  (3.9)  and  (3.10) 
directly  at  the  expense  of  arriving  at  a  pair  of  considerably  more 
complex  functional  equations  for  s^  and  t^ .  After  some  rather  tedious 
manipulation,  one  can  obtain  the  following: 


sin(a  -  2i|>)  j  [1  -  n  sin  a]  S^ct  -  ^)  -  S^Utt  -  a  -  \|i) 


in  a  |[1 


+  n  sin(a  -  2tp)]  S^(a  -  2tt  -  ip)  -  s^(-a  +  3if >)  >  (3.23) 


-  '--VvSI 


-49- 

[n  -  sin  a]t^(a  -  <f»)-s1na  C2ir  -  a  -  ip ) 

=  [n  +  sin(a  -  2<ji)]t^(a  -  2tt  -  +  sin(a- 2^)t^  (-a  +  3<|<)  .  (3.24) 

Similar  expressions  can  be  derived  for  s2  and  t2  ,  or  they  can  be 

obtained  directly  from  Sj  and  t1  ,  respectively. 

It  is  evident  that  (3.23)  and  (3.24)  are  linear  third-order 

functional  equations  of  the  unknowns  s  and  t  .  Theoretically, 

l  l 

solutions  to  this  type  of  equation  exist,  but  there  is  no  straight¬ 
forward  means  of  deriving  them,  except  in  certain  special  cases  [36]. 
This  is  analogous  to  the  task  of  finding  solutions  to  general 
differential  equations  of  order  greater  than  one. 

To  further  emphasize  the  complexity  of  equations  (3.23)  and 

(3.24) ,  consider  the  corresponding  functional  equation  developed 
by  Mai iuzhi nets  in  solving  the  scattering  by  an  impedance  wedge, 
namely 

(1  -  n  sin  a)s(a  -  <p)  =  (1  +  n  sin  a)s(a  -  2ir  +  t|i)  .  (3.25) 

Equation  (3.25)  is  a  linear  first-order  functional  equation  for  s(a). 
Two  important  differences  exist  between  (3.25)  and  either  (3.23)  or 

(3.24) : 
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1.  The  functional  equations  for  the  resistive  wedge  (3.23) 
and  (3.24)  are  of  third  order,  while  (3.25)  for  the  impedance  wedge 
is  of  first  order.  This  difference  in  itself  is  sufficient  to 
preclude  the  possibility  of  finding  a  solution. 

2.  Equation  (3.25)  is  a  special  type  of  functional  equation 
known  as  a  difference  equation,  where  the  functional  dependence  of 
the  argument  of  s  is  in  the  form  of  a  finite  increment.  In  (3.25), 
the  increment  is  2ir  -  2<i.  However,  (3.23)  and  (3.24)  are  not 
difference  equations.  This  is  evident  from  the  fact  that  the 
arguments  of  s,  and  t  contain  dependences  on  both  +a  and  -a, 
which  cannot  be  expressed  via  a  simple  increment.  Thus  the  author 
was  unable  to  apply  techniques  appropriate  to  difference  equations 
toward  finding  a  solution. 

In  order  to  possibly  obviate  the  difficulties  discussed  above, 
the  method  of  symmetric  and  antisymmetric  components  outlined  in 
Section  2.2  is  considered.  By  writing 


V 

1  I  ikpcosa 

"  2iri  e 

Sfi(a  -  $) 

.  . 

J  r 

te(a  -  <(>) 

(3.25) 


and  similarly  for  U°,  V°,  and  then  applying  the  boundary  conditions 
(2.24)  and  (2.25),  four  sets  of  functional  equations  are  developed: 


E-polarization  (symmetric) : 

n  sin  a£se(a  -  0)  -  se(o  -  0)3  -  2se(a  -  0) 

12  1 

=  -n  Sl'n  a[se(2ir  -  a  -  0)  -  se(2ir  -  a  -  0)]  -  2se(2ir  -  a  - 
1  2  1 

S?(a  -  0)  -  S6(a  -  0)  =  se(2ir  -  a  -  0)  -  se(2ir  -  a  - 

1  2  l  2 

se(a  -  ir)  *  -  se(ir  -  a) 

1  1 

Se(a)  =  -Se{2ff  -  a) 

2  2 

E-polarization  (antisymmetric): 

n  sin  «[s°(a  -  0)  -  s°(a  -  0)]  -  2s°(a  -  0) 

3  -n  sin  a[s°(2ir  -  a  -  0)  -  s°(2tt  -  a-0)]  -  2s°(2ir  -  a  -  1 
12  1 

S°(o  -  0)  -  s°(d  -  0)  »  s°(2tt  -  a  -  0)  -  s°(2*  - 

12  1  2 


S°(a  -  ir)  *  S°(ir  -  a) 
1  1 

S°(a)  ■  S°(2rr  -  a) 

2  2 


(3.27a) 

(3.27b) 

(3.27c) 

(3.27d) 

(3.28a) 

-  0) 

(3.28b) 

(3.28c) 


(3.28d) 
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H- polarization  (symmetric): 

2  sina  te(a  -  <|>)  -  n[te(a  -  ^ >  -  te(a  -  *)] 

1  1  2 

=  -2  sina  te(2ir  -  a  -  ip)  -  n[te(2u  -  a  -  <|i)  -  te(2ir  -  a 

1  12 

te(a  -  \p)  -  t^(a  -  <//)  =  —  t e( 2tt  -  a  -  ip)  +  te(2ir 

12  1  2 

te(a  -  ir)  =  -tfi(ir  -  a) 

1  1 

te(a)  =  -te(2ir  -  a) 

2  2 

H-polarization  (antisymmetric): 

2  Sin  a  t°(a  -  „,)  -  n[t°(a  -  *)  -  t°(a  -  4-)] 

1  1  2 

=  -2  sin  a  t°(2ir  -  a  -  <|i )  -  n[t°(2ir  -  a  -  ip)  -  t°(2ir  - 

1  1  2 

t°(a  -  ip)  -  t°(a  -  ip)  =  -t°(2ir  -  a  -<{')+  t°(2ir  - 

12  1  2 

t°(a  -  tt)  *  t°(ir  -  a) 

t°(.a)  3  t°(2ir  -  a) 

2  2 


*)]  (3.29a) 

a  -  ip) (3.29b) 

(3.29c) 

(3.29d) 

-  *)]  (3.30a) 

-  *)  (3.30b) 

(3.30c) 


(3.30d) 
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Again,  the  subscripts  1  and  2  refer  to  the  values  of  the 

corresponding  function  in  Regions  1  and  2.  As  noted  in  Section  2.2 

when  deriving  the  boundary  conditions,  the  only  difference  between 

(3.9)  and  (3.27-28)  or  (3.10)  and  (3.29-30)  is  in  the  final  two 

equations,  which  are  considerably  in  form.  However,  when 

6  0 

the  algebraic  manipulations  are  performed  which  eliminate  s  ,s  , 

2  2 

P  A 

t  ,t  as  was  done  for  equations  (3.23)  and  (3.24),  there  appears 

2  2 

to  be  little  gain.  The  resulting  equations  are: 

E-polarization: 


sin(a  -  24) 


=  -sin  a  | 


1 1  -  n  sin  a]  s^(a  -  4)  +  se(a  -  2ir  +  4) 


[1  +  n  sin(a  -  24))]  s^®  -  2tt  -  4)  +  S®(a  -  34)) 


sin(a  -  24;)  j  [1  -  n  sin  a]  s^(a  -  4)) 


-  s°(a  ■  +  4)  | 


=  -sin  a  |[1  +n  sin(a  - 


24)]  S°(a  -  2ir  -  4)  -  S  (a  -  34) 
l  1 


] 


(3.32) 


(3.33) 


H-polarization: 


[n  -  sin  a]  te(a  -  4)  +  s^n  ®te(a  -  2tt  +  4) 

1  1 

=  [n  +  sin(a  -  241)]  t^  (a  *  2ir  -  4)  -  sin(a  -  2\p)  t^(a  -  3^)  (3.34) 

[n  -  sin  a]  t°(a  -  4)  “  sin  a  t°(a  -  2ir  +4) 

1  1 

=  [n  +  sin(a  -  24)]  t°(a  -  2-ft  -  4)  +  sin(a  ■  24)t°(a  -  34)  (3.35) 

1  1 
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Similar  expressions  are  available  for  s^’0,!®’0,  or  they  can  be 

derived  directly  from  s®,0,t®’°  respectively.  Comparison  of 

(3.32-35)  with  (3.23-24)  indicates  that  the  decomposition  into 

symmetric  and  antisymmetric  components  has  done  little  to  simplify 

the  form  of  the  functional  equations.  There  has  been  no  reduction 

in  order;  (3.32-35)  are  still  third-order  equations.  If  any 

improvement  has  been  made,  it  is  that  (3.32-35)  are  now 

"difference"  equations,  in  that  the  arguments  of  the  unknowns  are 

of  the  form  a  +  An>  n  =  0,1, 2, 3,  where  An  is  a  finite  increment, 

and  the  increments  are  not  uniform  (i.e.,  A  .  -  A  ^  constant). 

n+i  n 

It  is  perhaps  in  doubt  whether  any  of  the  techniques  available 
for  solving  difference  equations  are  applicable  to  equations  of 
this  type. 

Unique  solutions  for  s®,0,t®’°  are  guaranteed  by  requiring 
that  they  satisfy  the  order  relation  (3.20)  with  a  >  0.  The 
regularity  condition  (3.22a)  must  be  modified  to  account  for  the 
presence  of  additional  plane  waves  in  the  incident  field  (see 
equation  (2.22));  the  result  is  that 


Se’°(a  -  *) 
1 

tf’°(a  -  it ) 


(a  -  2tt  +  $  )  *  +  (a  -  <f>  )  1 


must  be  regular  in  the  strip  ^  <_  Re  a  <_  2ir  -  <|» .  The  minus  (plus) 
sign  corresponds  to  the  symmetric  (antisymmetric)  components, 
respectively. 
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Equations  (3.32-33)  are  the  fundamental  functional  equations 
for  the  resistive  wedge  problem  formulated  via  Maliuzhinets' 
method.  The  author  has  attempted  to  solve  them  through  appropriate 
substitutions  and/or  factorizations,  with  no  success.  Methods 
described  in  texts  on  functional  equations  [36]  are  not  applicable 
to  third-order  equations  of  this  type.  Though  an  exact  solution 
has  not  been  found,  equations  (3.32-33)  are  amenable  to  iterative 
techniques  for  generating  approximate  solutions.  However,  it  is 
not  within  the  scope  of  this  work  to  pursue  those  methods  here. 
Instead,  an  alternative  formulation  based  upon  the  Kontorovich- 
Lebedev  transform  is  presented  in  the  next  chapter. 


CHAPTER  IV.  THE  KONTOROVICH-LEBEDEV  TRANSFORM 


4.1  The  Kontorovich-Lebedey  (K-L)  Transform  and  Its  Relationship 
to  the  Maliuzhinets  Representation 

In  1938,  two  Russian  authors,  M.  J.  Kontorovich  and  N.  N. 
Lebedev, put  forth  a  cylindrical  or  radial  transformation,  along  with 
the  corresponding  inverse,  which  is  analogous  to  the  LaPlace  (or 
Fourier)  transform  in  Cartesian  or  linear  coordinates  [21].  The 
transform,  which  now  bears  their  names,  found  applications  in  boundary 
value  and  diffraction  problems  where  the  unknown  functions  are  defined 
along  the  radial  coordinate  p  of  a  cylindrical  coordinate  system 

(p  » 2  )  . 

If  the  K-L  transform  of  a  function  f(p)  is  denoted  by 

f(v)  =  K[f(p)]  ,  (4.1) 

where  v  is  the  transform  variable,  then,  similarly,  the  inverse 
transform  is  written  as 


f(p)  =  K’^ftv)]  ,  (4.2) 

provided  both  K[f(p)]  and  K_1[f(v)]  are  defined. 

In  [21],  the  transformation  is  presented  in  the  form  of  a 
theorem  providing  sufficient  conditions  on  a  function  f(v)  for  the 
inversion  (4.2)  to  exist.  The  theorem  can  be  stated  as  follows: 
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Theorein  2: 

1 6 

Given  a  complex  number  k  =  |k|e  ,  0  <  6  <  tt,  and  w(v)  a  function 
of  the  complex  variable  v  *  a  +  it,  with  w(v)  satisfying  the  following 
condi tions: 

1.  w( v )  is  regular  (analytic)  in  the  strip  [ Re  v|  <  B»  6  >  0, 

2.  w(v)  is  an  even  function  of  v,  i.e.,  w(\>)  =  w(-v), 

3.  the  integral 


+  1x)w(a 


.  6t+(1t/2)(M-t) 

+  it  )  ( e 


dx  <  °° 


for  all  I  a |  <  6,  and 

4.  |<„  ♦  <<).{.  ♦  1t)|  eS,t(’/2)(l'l'T>  -  0  as  M  »  - 

uniformly  for  all  | o j  <  B. 

Then  supposing  that 


f(p) 


uw(y)e'i(ir/2)lJJ  (kp)dp 


where  p  >  0,  it  follows  that 


w(v) 


ei(ir/2)v 


H^}(kP) 


P 


for  all  v  in  the  strip  |  Re  v|  <  6. 

The  function  J  and  are  the  Bessel  function  and  Hankel 

v  v 

function  of  the  first  kind,  respectively. 

From  the  results  of  this  theorem,  it  is  convenient  to  define 


the  transformed  function 
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f(v)  « 

whereby  it  follows  that 

f(v)  -  K[f(p)] 

f(p)  -  K_1[~f(v)] 

i® 

v  sin  vir  f(v)e,vir  Hv‘y(kp)dv  ,  (4.4b) 

i°° 

the  second  integral  being  a  result  of  the  properties  of  Jv  and  H^1  ^ . 

It  is  important  to  point  out  that  the  transformation  put  forth 
in  the  manner  of  Theorem  2  must  be  used  with  caution.  By  assuming 
properties  for  f(v)  ( 1 . e . ,  w(v)),  and  then  defining  f(p)  via  (4.4b),  the 
theorem  restricts  the  class  of  functions  for  which  the  transform  (4.4a) 
exists.  In  practice,  it  is  the  properties  of  f(p)  which  are  known, 
and  therefore  it  is  necessary  to  show  that  f(v)  exists  as  defined  by 
(4.4a),  and  that  the  expression  on  the  right-hand  side  of  (4.4b)  does 
indeed  return  the  function  f(p). 

To  this  end,  Jones,  In  a  recent  paper  [22],  derives  sufficient 
conditions  for  f(p)  such  that  the  transform  Integral  and  its  inverse 
(4.4a,b)  exist.  In  particular,  Jones  defines  the  Inverse  as 


w(v)  e 


-i(ir/2)v 


(4.3) 


•/' 


*SL 


f (p  )  v  *  ( kp )  p 


(4.4a) 


vf(v) J^(kp)dv 


■1= 
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f(p) 


K_1[f(v)] 


eev  v  f(v)J^(kp)dv  (4.4c) 


in  order  to  guarantee  Its  existence  under  the  assumed  properties  of 
f(p).  An  outline  of  Jones'  theorem,  and  its  relation  to  the  properties 
of  the  functions  of  p  expected  for  the  resistive  wedge,  are  given  in 
Appendix  A.  Henceforth,  all  transformations  and  inverses  in  this 
chapter  will  be  assumed  to  exist,  either  in  the  sense  of  (4.4b)  or 
(4.4c),  based  upon  the  results  of  the  appendix. 

It  is  interesting  to  note  that  under  certain  conditions,  there 
exists  a  relationship  [24]  between  the  K-L  transform  of  f(p)  and  the 
Mai iuzhi nets  representation 


f(p) 


eikpcosa 


s(a 


-  <J> )  da 


(4.5) 


Y 

where  the  explicit  dependence  of  f(p)  on  has  not  been  shown,  and 
where  y  is  the  contour  of  Fig.  3.1.  Using  the  symmetry  of  the  contour 
Y  =  Yj  +  y2,  and  letting 


f (a)  *  j  [s(a  -  <f>)  -  s(2tr  -  a  -  $)]  ,  (4.6) 


(4.5)  becomes 


f(P) 


1 

STT 


eikpcosa 


f(a)da 


(4.7) 
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As  outlined  In  Appendix  A,  if  f (p )  =  0(p"1+a  e^p),  where  a,b  >  0,  then 
the  integral  equation  (4.7)  has  a  unique  solution  within  the  class  of 
functions 


f(«)  =  0(e(,-a)lI"al)  •, 


Im  a  I 


where  f(a)/sin  a  is  analytic  within  and  on  y,  given  by 


f(a) 


ik  sin 


f(p)  e“ikpC0Sa  dp  .  (4.8) 


Note  from  (4.6)  that  f(a),  as  well  as  y,  are  odd  about  a  =  ir,  i.e.. 


f(a)  «  -f(2ir  -  a)  , 


implying  from  (4.7)  that 


f(p) 


e1" k.pCOSa  “f(a)da 


Y 

1 


(4.9) 


Maliuzhinets  has  shown  [51]  that  in  the  event  the  constant  a  >_  1 , 
or  equivalently,  that  f(a)  is  bounded  at  the  end  points  of  y,  then 


f  (0)  =*  Zif(1») 


(4.10) 


4 


is  a  plane  wave  incident  upon  the  wedge  at  an  angle  $0  from  within 
Region  1  (see  Fig.  4.1).  The  included  angle  of  the  wedge  is  2i|». 
The  fields  (.11, V)  satisfy  the  Helmholtz  equation 


(v2  +  k2 ) 


U 

V 


3  0 


(4.14) 


in  Regions  1  and  2,  as  well  as  the  boundary  conditions  (2.6),  (2.7), 
the  radiation  condition  (2.9),  and  the  edge  condition  (2.12).  For 
the  time  being  is  assumed  that  Im  k  >  0.  Equation  (4.14)  is  also 
satisfied  by  the  scattered  fields  (US,VS). 

Throughout  the  remainder  of  this  analysis,  it  is  convenient  to 
define  the  modified  scattered  fields  (denoted  by  lower  case  letters) 


u(p  ,d> )  3  us(p  ,<j>)  +  ce*kp 

(4.15a) 

v(p  ,4> )  =  Vs(p,$)  +  de1lcp 

(4.15b) 

where 


»  * 

c 

us(o,*) 

d 

Vs(0,*) 

(4.16) 


The  quantities  c,d  are  independent  of  <j>  within  each  of  the 
Regions  1  and  2,  which  simply  means  that  the  scattered  field  is 
uniquely  defined  at  the  apex  of  the  wedge  within  each  region.  As 
was  discussed  in  Section  2.1  regarding  the  edge  condition  (2.12), 
it  is  possible  for  d  to  take  on  two  values,  one  each  in  Regions  1 


Region  1 


Region  2 


\ 

71  \ 


(j>=2ir_i|/ 


Fig.  4.1:  Geometry  for  the  Application  of  the  Kontorovich- 
Lebedev  Transform  to  the  Scattering  by  a  Resistive 
Wedge. 
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and  2,  but  due  to  the  continuity  of  l)  across  the  boundary,  the  value 
of  c  is  the  same  in  both  regions.  From  equation  (4.16),  it  is 
evident  that  (u,v)  are  zero  at  p  =  0. 

Expressing  (4.14)  in  cylindrical  coordinates,  and  substituting 
(4.13)  and  (4.14),  the  following  equations  for  (u,v)  are  obtained: 


2  32  + 

Pz  -  + 

3P2 


3p 


—  +  k2p2 

'  U  ' 

=  ike  e,k(> 

■c  * 

3<j>2 

V 

d 

Assuming  the  additional  constraints 


11m’  fe 

p  o 


u 

V 


L  vj 


(4.17) 


(4.18) 


which  are  consistent  with  the  edge  and  radiation  conditions  for  p  and 
^-directed  components  of  the  scattered  fields,  then  application  of  the 
K-L  transform  to  (4.17)  leads  to  a  differential  equation  for  the 
transforms  (u,v)  [21]; 


di_  +  v2 

u  (v  ,4> ) 

c 

d4>2 

v(v,<t>) 

sin  vtt 

.  d  . 

(4.19) 


The  general  solution  to  this  equation  can  be  written  in  the 

form 


U(\M$) 

A(v) 

COSvij)  + 

B(v) " 

-In  ,  2<  n'iU/2)v 

’  c  ” 

v(v  »(^  ) 

•  « 

C(v) 

•  a 

D(v) 

Sln  +  V  Sin  VTT 

d  _ 

(4.20) 


-65- 


where  the  arbitrary  functions  A,B,C,  and  D  are  independent  of  *.  The 
task  is  to  find  specific  values  of  A,B,C,  and  D  such  that  (u,v) 
satisfy  the  appropriate  transformed  boundary  conditions.  In  terms 
of  (u,v),  the  conditions  are: 


E-Polarization: 


3Ut 
- 1 

2ikp 

+ - *=•  u 

3<j> 

3(j) 

n  1 

2iko  j-eikpcos(ip-(j>0)  _  ceikP] 


4,=^, 


(4.21a) 

(4.21b) 


♦ 


34> 


3U 


=2ir-ip  ^ 


(|)  =— Ip 


-  21k.u 
n  1 


^*2ir-iJ) 


=  e~i kpcos (ip-Hj>0) _  ceikp  J 


(4.21c) 


U  “  u  I 

1  <j>  =2tt— ip  2  1 4>  =— ip 


0  . 


( 4 . 21 d ) 
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H-Polarization: 


3<j> 


iknp 

O 

v. 

1 

=  -ikp 


sin  (*  -  *o)e-ik,,“s(*-'*o) 


-  J  (d  -  d  )e 

4  1  2 


ikp 


(4.22a) 


3  V. 


34> 


4>=>l> 


3V, 

3<)> 


< P ** 


(4.22b) 


<f>=2ir-i|; 


ikrip 

2 


v 

l 


t(/=2n-iji 


v 

2 


4>=-'i' 


=  ikp 


sin(ip  +  * 


jgikpcos(ip+  <t>0) 


-i'V 


d  )e 

2 


ikp 


(4.22c) 


3  V 


1 


3  V 


2 


3<j> 


^=2ir-i{) 


0 


(4.22d) 


The  K-L  transform  of  equations  (4.21)  and  (4.22)  are  obtained 
directly  by  multiplying  each  equation  by  (l/p)H^  (kp)  and  integrating 
along  (Q,<»).  The  results  are: 
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E-Polarization: 


du 

1 


_  d%_ 

+  2i_k 

1 

I  u 

d<** 

n  j 

1  il 

0 

41  e'i<,,/2)v 

Pm 

n 

sin  vir 

u 

_ 

l 

H^OciJdp 


sin  \)(tt  -♦  +’!») 
sin($0  -  *} 


=  0 


<t>='l» 


du 

HfLi 


<ps2ir-<p 


du 

•  as2" 


-i p 


-  ^  f  u  H(1)(kp)dp 

n  Jn  1  2*-*  V 


4i  f  sin  V(1T  -  4>0  -  ’!») 

~  sin  vir  CNJ  "  sin(<(>0  +  4>) 


-  u 

$=2ir-i|i  2 


*  0 


♦  “■'I' 


H-Polarization: 


dv. 


<t>=^ 


up-  f  (v  |  -  v  H" 


(kp)dp 


-i(tt/2)\J  r 

21  \in  V7T  [  J  {di  -  d2)v  +  sin  v{lT  •  *o  +  *} . 


4>=V< 


(4.23a) 

(4.23b) 

(4.23c) 

(4.23d) 

(4.24a) 


0 


(4.24b) 
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dv  dv 

_ _  _ 2 

^  (j>=2ir-ijj  ^  $=-41 


(4.24d) 


Some  note  regarding  the  regions  in  the  complex  v  plane  for  which 
equation  (4.20),  and  boundary  conditions  (4.23)  and  (4.24)  are  defin 
is  necessary  at  this  point. 

For  (u,v)  =  0(p  ),  6  >  0,  it  is  shown  in  Appendix  A  that  (u,v), 
and  hence  the  left-hand  side  of  (4.20),  are  defined  for  [Re  v|  <  5. 
Inasmuch  as  the  unknowns  A,B,C  and  D  are  determined,  the  RHS  of  (4.20) 
may  or  may  not  provide  an  analytic  continuation  of  (u,v)  into  any 
additional  regions  where  it  is  properly  defined. 

From  similar  considerations,  the  LHS  of  (4.23)  and  (4.24)  are 
also  defined  for  [Re  vj  <  6,  while  the  RHS  are  analytic  in  ] Re  v|  <  1, 
having  poles  at  v  equal  to  a  nonzero  integer. 

It  is  important  to  note  the  various  regions  of  analyticity  of 
the  transformed  quantities,  since  the  inversion  contour  of  (4.4b,c) 
must  lie  within  the  intersection  of  these  regions. 

With  these  caveats  in  mind,  equation  (4.20)  may  be  substituted 
into  the  transformed  boundary  conditions  in  an  attempt  to  determine 
the  unknown  coefficients  A,B,C,  and  D  (each  defined  in  Regions  1  and  2). 
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From  (4.23b,d), 


A  (v)  cosviii  +  B  (v)  sin  -  A  (v)  cos  -  B  (v)  sin  vxp  =  0 
l  l  2  2 


implying 


fc= 

CM 

-  ^)  +  B  (v)  : 

sin 

v(2ir  - 

*) 

-  A 

2 

(v)  COS 

Vlp 

+  B 

2 

(v) 

sin 

=  o  , 

A  (v) 

cos  v(tt  - 

Jd 

[A  (v) 

COS 

\>*TT 

+  B  (v) 

sin  vtt ] 

(4.25a) 

2 

COS  Vtp 

l 

l 

B  (v) 

2 

sin  v (tt  — 
sin  \>\p 

Jd 

[A  (v) 
l 

sin 

VTT 

-  B  (v) 

1 

COS  vir] 

(4.25b) 

Similar  results  can  be  derived  from  (4.24b,d),  i.e., 


Mv)  =  -  51  Vfn^"^  cos  V1T  +  Vv)  sin  vir^  (4-26a) 

Mv)  =  '  ^COS^'^  CC!(V)  Si"  V1T  ’  °1(V)  C0$  V1T]  (4>26b) 


Substitution  of  (4.25,  26)  into  the  remaining  boundary  conditions 
(4.23a,c)  and  (4.24a ,c)  leads  to  the  four  rather  complicated  equations 
below: 

E-Polarization: 

cos  visin' v<i>  [-sin  vU  -  *}Vv)  +  cos  v(lT  -  +} 


2  i  k 


<t>=^ 


H^(kp)dp  =  K(  v  ,ip ) 


(4.27a) 
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v  sin  vir 
cos  \h|j  sin  vifi 


[-sin  v(tt  +  <ji)  A^(v)  +  cos  v(it  +  4/)  B^v)] 


2i  k 

n 


$=2it-<|i 


H^(kp)dp  3  -K(v,-i|0  (4.27 


H-Polarization: 


-v[C  (v)  sin  viJj-  D^(v)  cos  vi|j]  +  -^!p 


v  -  v 

1 1  ip  ^ 


•  H^(kp)dp  =  L(^,'k)  (4.2f 


-v 


(v)  sin  v(2ir  -  <|>)  -  0  (v)  cos  v(2tt  .  4,)] 


ikn 

•  ~r 


4>  =2tt 


L(v, 


(4.2 


The  functions  K(v,i|i)  and  L(v,<| >)  are  given  by 


K(v  ,i|)) 

4i_ 

e-i  U/2)v 

r\\  mm  . 

n 

sin  vtt 

wv 

• 

L(v,i)<) 

=  2i 

e-i  (it/2)v 

sin  vir 

sin  v(ir  -  4>0  +  4>) 
sinUQ  -  *J 


.  (4.2 


5-  (d  -  d  )v  +  sgn  t|»  sin  vU  -  <j) 
2  i  2  1 


(4. 


Equations  (4.27)  and  (4.28),  defined  in  the  strip  (Re  v|  < 
min  (5,1),  are  a  pair  of  coupled  equations  for  the  remaining  unknowns 
A  ,B^ ,C^ ,  and  ,  and  their  solution  provides  all  that  is  necessary 
for  determining  the  unknown  transformed  fields  (u,v).  However,  these 
equations  cannot  be  solved  explicitly  in  their  present  form,  due  to 
the  presence  of  the  "untransformed"  fields  (u,v)  under  the  integral 
signs. 

Two  methods  for  simplifying  (4.27)  and  (4.28)  by  eliminating 
the  "untrans formed"  fields  are  presented  in  the  following  subsections. 

In  4.2.1,  an  approach  leading  to  second  order  difference  equations  is 
developed.  The  technique  parallels  that  used  by  Lebedev  and  Skal'skaya 
[23]  in  solving  the  impedance  wedge  problem.  A  second,  more  general 
method  leading  to  Fredholm  integral  equations  of  the  second  kind  for  the 
unknowns  is  described  in  4.2.2.  The  method,  developed  by  the  author, 
allows  the  well-established  theory  of  linear  integral  equations  to  be 
applied  toward  solving  the  equations. 

4.2.1  Difference  Equation  Method 

The  method  to  be  presented  requires  that  one  of  the  following 
two  conditions  is  satisfied  by  the  unknown  functions  (u,v); 

0 )  5  >  1 ,  or 

(2)  (u,v)  may  be  analytically  continued  into  the  region 
I  Re  v|  <«',«'  >  1. 

The  net  result  of  either  of  these  requirements  is  that  (u,v)  or  their 
continuations  are  analytic  in  the  strip  |Re  vj  <  1  +  e>  e  >  0.  For 
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the  moment  it  will  be  assumed  that  such  is  the  case.  Whether  it  is 
indeed  true  will  be  discussed  in  another  part  of  this  work. 

Equations  (4.27)  and  (4.28)  are  now  multiplied  by  e^ir>^v, 
and  then  v  is  replaced  first  by  (1  +  v)  and  then  by  (1  -  v).  The 
resulting  set  of  equations  in  (1  +  v)  have  a  common  region  of 
analyticity  with  the  corresponding  set  of  equations  in  (1  -  v),  which 
from  the  above  requirements  is  given  by  | Re  v|  <  e.  By  subtracting 
one  set  from  the  other,  and  making  use  of  the  identity 


e1(./2)(l«v)  H(l)(ko)  .  et(»/2)(l-v)  „(l)(ko)  .  |ia.ei{»/2)v  „(> )(ko) 

(4.30) 

there  results  a  new  system  of  equations  for  the  unknowns  A* ,  Bj , 

Ci ,  Dx : 


E-Polarization: 


P+(l  +  v)F(l  +  v)  -  p+(l  -  v ) F(1  -  v)  +  q+(l  +  v)G(l  +  v) 

4  8i  COS  v (tt-4>  +if>) 

-  q+0  -  v)G{l  -  v)  -  -  COS  ip  F(v)  +  sin  vij)G(v)  =  —  - iTfTviT - 


(4.31a) 

p_(l  +  v)F(l  +  v)  -  p_(l  -  v)F(l  -  v)  +  q.(l  +  v)G(l  +  v) 

-  q  (1  -  v)G(l  -  v)  +  —  (cos  v(2it  -  if))F(v)  +  sin  \>(2n  -  <J>)G(\>)3 

“  n 

COS  v(ir  -  4>0  -  <J>) 
sin  vir 


-  8i 
n 


(4.31b) 


where 


P±tv) 

sin  vir  sin  v(ir  +  ip) 
sin  vip  cos  vip 

q+Cv) 

sin  vtr  COS  v(tt  +  ip) 
sin  vtp  cos  vip 

F(v) 

=  V  eiM)\(v) 

G  (v) 

=  v  ei(lT/2)vB  (v)  . 

H-Polarization: 


sin(l  +  v )4» F '  ( 1  +  \j)  -  sin(l  -  v)\pF'(1  -  v)  -  cos(l  +  v)tpG'(l 


+  cos ( 1  -  v)^G'{1  -  v)  +  r)[q+(v)F' (v)  -  p+(v)G'(v)] 

n(^0  ”  'l'^C0S  "  <J>0  +  ^ 

41  sin  vtt 


sin[(l  +  v)(2tt  -  ip)]F'(l  +  v)  -  si n[(  1  -  v)(2ir  -  ip)]F'(l  -  v) 

-  cos[(l  +  v)(2ir  -  ip)]G'(1  +  v)  +  cos [ ( 1  -  v)(2ir  -  if>)]G'(l  - 

sin4_  +  tp)cos  v(ir  -  <f>  - 

-  n[q_(v)F'(v)  -  p_(v)G'(v)]  3  41 - - — 

where 

F'(v)  »  e1(ir/2)v  C^v) 

G’(v)  *  gi (ir/2)v  D  (v)  > 

l 


+  v) 

(4.32a) 

v) 

^  (4.32b) 
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Equations  (4.31)  and  (4.32)  each  constitute  a  pair  of  coupled 
second  order  linear  functional  equations  for  F(v),G(v)  and  F'(v), 
G'(v),  respectively.  Although  one  may  observe  several  symmetries 
associated  with  the  equations,  there  is  no  straightforward  means 
for  finding  their  solutions  known  to  the  author. 

Certainly  the  task  of  developing  expressions  for  F  and  G 
(and  likewise  F'  and  G')  would  be  easier  if  the  functional  equations 
could  be  uncoupled. 

To  this  end,  the  decomposition  derived  in  Section  2.2  is 
applied  to  the  transformed  functions  (u,v).  From  the  equivalent 
of  equation  (2.18),  one  can  define 


•  ^  « 
u  (v,4>) 

ue(v  ,<(> )  +  u°(v,<fr) 

1 

1 

1  i 

V  (v,<f>) 

2 

ve(v  ,4> )  +  v°(v,*) 

l  1  J 

1  l  J 

where 


e,o 

u 

u  *  u  (v,2ir  -  <{>) 

l 

l  l 

3 

V  (v,4>)  ±  V  (v ,2ir  -  <j>) 

l 

L  i  i 

(4.33) 


(4.34) 


The  plus  (minus)  signs  correspond  to  the  symmetric  (antisymmetric) 
components.  It  is  convenient  to  write,  analogous  to  (4.20)  . 
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r 

e,o 

u 

l 

A*,0(v)" 

V’0(v)  ' 

2ie~i  (ir/2)v 

r  ce‘°  i 

= 

cosv<j>  + 

sin  \>$  +  - 

y  v  sin  w 

y6»0 

ce*°Cv) 

De*°(v) 

de*° 

.  1 

L  l 

L  i  J 

.  i 

(4.35) 


Substituting  (4.35)  into  (4.33)  leads  to  the  following  results: 


[  [Bj  +  B°]  . 


(4.36a) 

(4.36b) 


and  similarly  for  and  In  addition,  from  (4.34) 


ce  =  2c 
l 

de  *  2d 
i  l 

c°  =  d° 


From  either  (4.34)  or  the  transforms  of  ( 2 . 24c ,c ’ )  and  (2.25c,c'),  it 
is  easily  shown  that 


du; 

d<t> 


dv 
d*1 ! 


4>*tt 


0  , 


$“ir 


-  -i 
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F  (1  +  \>)  F  (1  -  j)  .  e.  ,  « 

- '  -  .  Q  ■  ■  ■  +  A  sin  -  Al  F  (v) 

sin  (1  +  v)^  sin  (1  -  v)\ji  n  sin  vir  ov  1 


16i  sin  "  4>0)si n 
n  sin  vtt 


(4.39b) 


H-Polarization: 


sin[(l  +  vHtt  -  »)] 
sin(l  +  v)ir 


f;(i  +  v) 


sin[(1  -  vHi r  -  »)] 
sin(l  -  v)ir 


f;o  -  v> 


f;<v) 
sin  \>ii> 


iTn^v?  Csin(^o  ‘  *>cos  v(ir  -  *0  +  *> 


-  sin  (<(>o  +  i|i)cos  v(rr-  $o  -  i(;)]  (4.40a) 


ces[(l  +  v)(n  -  »)]  F,(1  +  x  .  cos [( 1  -  v)U  -  »)]  F,n  .  % 

sin(l  +  4ir  V  v;  si n(  1  -  vjn  V  ; 


+  n 


fjM  _  _ 

cos  vij/  sin  vtt 


4i 


[si n(<|>  -  4* )cos  v (it  -  $  +  i|») 


+  sin(4>0  +  4«)cos  v(n  -  <f>0  -  41)]  (4.40b) 


Equations  (4.39)  and  (4.40)  are  the  fundamental  functional 
equations  for  the  resistive  wedge  scattering  problem.  While  an 
equation  of  similar  type  has  been  solved  by  Lebedev  and  Skal'skaya 
in  [23],  the  author  has  not  been  able  to  generate  a  solution  using 
their  method. 
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A  general  solution  to  any  of  equations  (4.39,40)  can  be 
written  in  the  form 


F(v)  =  P  (v)F(l)(v)  +  P  (v)F(2)(v)  +  F  (v)  ,  (4.41) 

1  2  P 


where  F^,F^  are  the  two  independent  solutions  of  the  corresponding 
homogeneous  (RHS  =  0)  equation,  and  Fp(v)  is  a  particular  solution 

of  the  inhomogeneous  equation.  P  and  P  are  arbitrary  even  periodic 

1  2 

functions  with  period  unity. 

Determination  of  a  unique  solution  to  (4.41)  requires 
restricting  the  solution  to  a  particular  class  of  functions 
satisfying  a  specified  analyticity  condition,  along  with  a 
prescribed  behavior  as  |I m  v|  -*•  <*>. 

More  specifically,  it  is  possible  to  write 


F  (v)  *  -4ic  +  ve  ^  sin  vtt 


=TT 


H^dcp)  jj2-  (4.42a) 


FQ(v)  *  ve  sin  vtt 


<()=7r 


H(l)(kp)  &- 
v  p 


(4.42b) 


with  identical  expressions  for  F‘  in  terms  of  v  ,  d  ,  and  for  F'(v) 

e  ii  o 

in  terms  of  v°,  where 

l 


r  e 

0  i 

u 

=  1 
2 

u 

+ 

u 

e 

0 

.  V  . 

L  v 

+ 

V  _ 

(4.43) 


1 
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From  (4.39-4.43)  and  behavior  of  (u  .  v  |  )  at  p  =  0 

i  <P  — ir 

and  for  p  -*■  «,  the  following  conditions  are  satisifed  by  Fg,F  ,Fg,F^; 

1.  Fe,F^  are  even  functions  of  v,  while 

F1,F„  are  odd  functions  of  v. 
e  o 

2.  The  function  Ffi(v)  is  analytic  in  the  strip  (Re  v|  <6, 

6  (5 

where  u  =  0(p  ),  5  >  0,  for  p  •+■  0.  Similar  results  hold  for 

l 

F  ,F',F'  in  terms  of  u°,ve,v°,  respectively, 
o  e  o  ill 

3.  The  function  Ffi(v)  approaches  zero  when  | Re  v |  <6  and 

llm  v|  -*■  «,  provided  ue|.  contains  no  geometrical  optics  scattered 

^  <J>  =1T 

fields  (recall  the  decomposition  (2.8)),  which  is  always  true  for 
4>0  -  2\p  >  0-  Similar  results  hold  for  FQ,Fg,  and  F^.  In  [23], 
Lebedev  and  Skal'skaya  write  Fe(v)  =  F®9(v)  +  Fed(v),  equivalent 
to  (2.8),  and  then  state  that  Fed(v)  -*•  0  for  | Re  v|  <  5, 

|Im  v|  -*■  »  for  all  This  *s  unnecessary,  since  by  solving  the 

problem  for  <t>Q  -  2<|i  >  0,  (whereby  F69  *  0),  the  solution  for 
♦  -  2*  <  0  can  be  obtained  via  an  appropriate  analytic  continuation. 
Further  discussions  regarding  the  behavior  as  | Im  v|  »  can  be  found 
in  Appendix  A. 

Conditions  (1-3)  above  are  sufficient  to  uniquely  determine 
the  solution  (4.41 ) . 

When n  =  0,  equations  (4.39)  can  be  solved  directly,  yielding 
the  known  solutions  for  a  perfectly  conducting  wedge  (E-polarization) 


F«(v) 


nao 


cos  v(ir  -  ) cos  vq» 

-4i  - t — - — « - 

COS  v(tt  -  ip) 


F0(») 
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nao 


sin  v(7T  -  $0)s1n  v\p 
sin  v(it  -  ip) 


(4.44a) 


(4.44b) 


Although  the  solution  to  (4.40)  (H-polarization)  with  n  =  0 
cannot  be  as  easily  deduced,  the  resulting  equations  can  be  solved 
via  Fourier  transform  techniques,  yielding  the  correct  expressions. 


Furthermore,  it  is  relatively  straightforward  to  solve  (4.39- 
40)  for  the  special  cases  ip  =  0,  ir/2,  corresponding  to  a  resistive 
half-plane  and  full  plane,  respectively.  Both  solutions  are  known, 
the  latter  being  the  simple  geometrical  optics  reflected  and 
transmitted  plane  waves,  while  the  former  having  been  given  in 
[27,28]. 

The  lack  of  available  methods  for  solving  functional  equations 
of  the  type  (4.39-40)  has  prevented  the  author  from  obtaining  an 
exact  solution  for  arbitrary  values  of  \p.  However,  since  (4.39-40) 
were  derived  from  (4.27-28)  based  on  the  assumption  that  (u,v)  were 
analytic  in  a  strip  |  Re  v|  <  <5,  with  <5^1,  it  is  reasonable  to  ask 
if  a  different  set  of  equations,  replacing  (4.39-40),  can  be  derived 
in  the  more  general  situation  where  6  >  0.  Such  a  derivation  is 
presented  in  the  next  section. 

4.2.2  Integral  Equation  Method 

The  starting  point  for  this  section  will  be  equations 
(4.27-28).  It  is  beneficial  to  immediately  express  them  in  terms 
of  symmetric  and  antisymmetric  components.  This  is  achieved  by 
substitution  of  (4.36-38),  then  adding  and  subtracting  (4.27a)  and 
(4.27b),  and  similarly,  (4.28a)  and  (4.28b).  The  net  results  are 
the  following: 


E-Polarization: 


^±tmetW2)»f‘  „o(kp)dc 

cos  vi|)  n  J  1  v 


=  CSC  VTT 


sin  v (tt  -  $_  +  t|i)  sin  v(ir  -  <p.  -  <p) 

pr  _ _  o  _  0 

sin(ir  -  4>0  +  t)  sin(ir  -  <f>0  -  it>) 


(4.45a) 


Fo^  2  i  k  i  (if/2)v 

sin  \>\p  ~  ri  e 


k.X1,(kp,d° 


4i 


CSC  vir 


sin  v(ir  -  4>  +  iji)  sin  v (tt  -  4>0  -  ) 


sin(:r  -  <t>0  +  tl>)  sin(ir  -  4>Q  ~  <p) 


(4.45b) 


H-Polarization: 


sin  v(ir  -  »)  F, .  s  jkn  i(n/2)v 

sin  vtt  V  '  2  e 


V  -  v|] 


.  .  H(l)(kp)dp 
>  =lj>  v 


=  2i  esc  viT[n(di  -  d^)  +  sin  v(tt  -  ij>0  +  iJj)  -  sin_v(ir  -  <t>0  -  ^)] 


(4.46a) 


cos  v(ir  -  »)  . 

sin  vit  ho'v' 


ikn  1(ir/2)v 
2  e 


/*  00 

J 

0 


♦  H^1 >(kp)dp 


=  -  2i  esc  vir[sin  v(ir  -  4>0  +  (|»)  +  sin  v(ir  -  $  -  ip)]  (4.46b) 


From  this  point  on  it  will  be  assumed  that  (j>Q  -  2ip  >  0 
(implying  $  <  ir/2),  which  for  reasons  stated  earlier  allows  Fg(v)  -*•  0 
for  |Im  v|  -*■  »,  | Re  v|  <  5.  Similar  behavior  is  exhibited  by 
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0 

In  [23]  it  was  shown  that  by  decomposing  u  in  geometrical 
optics  and  diffracted  components,  i.e., 


u 


e 

l 


ueg  +  ued 

l  l 


(4.47) 


the  transformed  quantity  u^  -  u®g  satisfies  all  the  conditions  of 
Theorem  2,  and  hence  the  following  representation  is  valid 


/> 


e  eg 

u  -  u  = 
l  l 


T_ 

4 


J  _io 


y  sin  yir  [u“(y) 


-  ue9(y)]eWl)(kp)dp  . 


(4.48) 


Therefore 


u^jH^fkpJdp  •  -T,\  H^l)(kp)- 

J  o 

sin  yir[ue  -  ueg]el7riJ  H^(kp)  dy  dp 


1  1 


(4.49) 


From  the  conditions  of  Theorem  2  along  with  results  from  [21],  it  is 
allowable  to  exchange  the  order  of  integration  in  (4.49).  With 
the  aid  of  the  identity 


ei(*/2)x 


h'1 >(k0) 


iT  K*Hke)  ' 


the  integration  on  p  can  be  expressed  in  terms  of  the  integral 
[52] 


(4.52) 
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A  00 

■i(ir/2)v  j  Kv C- i kp ) ( - i kp ) dp  =  e  ^T^^v[cos  pir  +  cos  vir]  1 


Re  ji  |  +  [  Re  v  j  <  1  ;  Ira  k  >  0  . 


The  net  result  Is 


r 

o 


[u*  -  ueg]  H^1 ^ (kp )dp 


1  e-i(if/2)v  • 
"  2k  6 


psin  pir[u®  - 

COS  prr  +  COS  vtt 


dp 


(4.50) 


Jones  has  shown  [22]  that  the  integral 


ug6(p) 


ugeH(l)(kP)  ^ 
1  ^  ^ 


(4.51) 


is  uniformly  convergent  for  Re  p  =0,  and  hence  (4.51)  may  be 
substituted  into  (4.50)  and  the  order  of  integration  exchanged; 
therefore 


f 

*  n 


[u*  -  uSe]H<’>CKp)dp  - 


p  sin  pit  uxe 


r,e  J(w/2)p 


COS  PIT  +  COS  vir 


dp 


JL  p"*  U/2)v 
2k  e 


jge 

i  p 


r 

<k 

*'-i« 


v  j 
•1° 


...  =  WL(l)(»)i,hffl»  du 


COS  y7T  +  COS  vir  p 


(4.52) 
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By  means  of  the  identity 


u  e1(*/2)vH(i)(kp)  =  [ic^f-ikp)  -  Ku+1C-ikp)]  , 


it  follows  that 


loo 

/ 


sin  yir  e1^2^  (l),  . 

cos  yrr  +  cos  wir  pHp 


-loo 


loo 

¥■  J 


-i« 


sin  pit _  _  . 

cos  hit  +  cos  vir  tKy_1(-ikp)  -  K^+i (-ikp)]dp 


k£ 

IT 


ioo-1 


-i®+i 


{  '  J 

-i<*>-i  i“+i 


—  sin,.M - K  (-ikp)dy  (4.53) 

COS  Utt  -  COS  Vir  yv 


The  function  K^(-ikp)  is  analytic  within  the  strip  |lm  y|  <1,  for 
Im  k  >  0,  and  |K^(-ikp)|  +  0  as  |  Im  y|  -*•  ®  within  the  strip. 
Therefore,  (4.53)  may  be  evaluated  using  residue  theory,  whereby  the 
terms  on  the  LHS  and  RHS  of  (4.52)  containing  u9e  cancel,  leaving 


ueH(I>(kp)dp  =  -l-e-1!*/2)v 
1  v  lK 


1o 


p  sin  Fir  ei(7r/2)tJ  ~ue{v)du 


COS  yir  +  COS  vir  i 

-i"  (4.54) 


|  Re  v  |  +  |  Re  y  |  <  1  ;  Im  k  >  0  . 


It  is  easily  shown  that 


-a  -1(ir/2)y 

U!(w)  "  y  sin  yV •  £Fe(w)cos  >*(*  -  ♦)  +  4ic^  » 


(4.55) 
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and  In  addition. 


/>  i® 


41c 


-i« 


COS  11TT  +  COS  VTT 


dy  = 


8cv 


Sin  VTT 


(4.56) 


Substitution  of  (4.54-56)  into  (4.45a)  results  in  an  Integral 
equation  for  Ffi(v),  viz 

/i® 

c,Pi  L(y  »),J?ps  vjp  (u)d 

cos  yir  +  cos  \>ir  M 


-1 


4i  cos  vij) 
n  sin  vir 


f  Sin  v(ir  -  4>0  +  <p)  Sin  v(ir  -  <pQ  -  <f-) 
s1n(<(>_  -  +  sin((ji  +  ?j 


(4.57) 


|Re  v|  +  | Re  y|  <  1 


Equations  analogous  to  (4.54-56)  exist  for  the  functions  u°,v°,ve,v0, 
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and  v°;  hence  it  is  possible  to  generate  similar  integral  equations 
from  (4.45b)  and  (4.46)  for  FQ,Fg,  and  F^. 

Prior  to  presenting  the  integral  equations,  it  is  convenient 
for  later  analyses  to  put  v  =  ii,  y  *  it1,  and  define 


[fgCx ) ,fQ(T ) D 


Cf;(T),f;(T)] 


[Fe(iT),F0(iT)] 


CF;(1t),F'(1t)] 


(4.58) 


In  terms  of  the  new  unknowns,  the  integral  equations  can  be  written 
as* 


E-Polarization: 


Ch (,7T  -  ^  >T  *  ch  lt»T 
ch  ITT1  +  ch  ITT 


fg(T 1  )dt 1 


4i  ch  jrr 
n  sh  itt 


sh(.ir  -  4>0  +  ip) t  sh(ir  -  <J>Q  -  \{i)t 
sin($.  -  *J  +  sin(4>  +  J] 


V*> 


s hh,  -  t)j.'  il?.  f  ,  ,)d  , 

Ch  TTT1  +  Ch  ITT  fou  )ar 


4^  sh  ijrr 
n  sh  irx 


sh(ir  -  4>0  +  t<#)t  sh(ir  -  i(>0  -  )x 
sin(<{»0  -  ip)  sin(<j>Q  +  *) 


H-Polarization: 


f;(x)  + 


0 

•f 


Sh  TIT1  Sh  TTT 


sMrp  sh(ir  -  i(|)t  ch  TTT1  +  ch  TTT 


dx ' 


*  4i 


c h ( ir  -  4>0)x  sh  ipr 


sh(ir  -  iJj)t 


Sh  TTT  '  sh  ITT 


ch  tJtt'  ch(.ir  -  ch  tt'  +  ch  ttx 


dx' 


Sh(ir  -  <f>Q)x  ch  tJtt 
Ch(ir  -  t(/)x 


The  functions  sh(  ),  ch(.  },  represent  the  hyperbolic  sine  and 
functions,  respectively. 


(4.59a) 


(4.59b) 


(4.60a) 


(4.60b) 

cosine 
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The  even  or  odd  property  of  each  unknown  has  been  used  to 
reduce  the  Integration  Interval  to  the  serai -Infinite  line  [0,«). 

Equations  (4.59-60)  are  the  fundamental  equations  for  the 
resistive  wedge  scattering  problem  formulated  via  the  K-L  transform 
technique.  The  author  is  unaware  of  any  derivation  of  this  type  In 
the  literature  on  the  transformation.  The  equations  are  Fredholm 
integral  equations  of  the  second  kind  (when  n  t  0). 

All  four  equations  can  be  written  in  the  form 


■{ 


K(t  »t  1 )  f  (t  '  )dx 1  = 


0  <  T  < 


(4.61) 


The  functions  f(x),  g(x)  are  simultaneously  either  even  or  odd 
functions,  while  K(x.x')  is  simultaneously  either  even  or  odd  in  both 
x  and  t'.  The  parameter  x  is  proportional  to  either  n  (H-polarization) 
or  rf1  (E-polarlzation) .  The  kernels  of  the  Fredholm  integrals  are 
all  bounded  in  the  quarter-plane  0  <_  x  <  «,  0  <_  x'  <  ».  In  addition, 
the  following  order  relations  hold 


|K(x,x')|  *  o(e"'*'l Re  Tl )  ;  | Re  x |  -*•  »  (4.62a) 


for x  ’  “constant.  Likewise 

|K(x,x’)|  -  0(e“^ir"'*'^Re  T*  ^ )  ;  | Re  x '  |  -*■  »  (4.62b) 


for  x  *  constant.  As  a  result,  it  can  be  shown  that 
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f  |  K(x  ,t ' )  |  dx  <  «  ;  for  real  t' 


rm 


| kCt ,T ' )  I  dr 1  <  «  ;  for  real  r 


(4.63) 


o 

i.e.,  the  kernel  K(t,t')  Is  absolutely  integrable  in  both  t  and  t'. 
However,  it  is  not  integrable  in  the  quarter-plane,  that  is  to  say 


rm  rm 

| K(t  ,r 1 ) | 2  dr  dr 1 

J  « 

o  o 

Is  unbounded.  Furthermore,  the  inhomogeneous  terms  g(r)  of  the 
integral  equations  (4.59-60)  satisfy 


| g(r) |  =  0(e’(V2*)|Re  x|  )  ;  |Re  t|  -  -  .  (4.64) 


Since  it  has  been  assumed  that  $  -  2<|>  >  0  for  this  analysis,  it  is 

also  true  that 


! 


00 


jg(r)|dr  <  »  . 


(4.65) 


o 

The  analytlcity  properties  of  the  kernels  and  Inhomogeneous  terms, 
represented  by  K(r,r')  and  g(r)  respectively,  are  also  easily  determined, 
and  can  be  summed  up  as  follows: 


1 


I 
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lemma  1 : 

There  exists  a  <$  >  0  such  that  g(t)  is  analytic  in  the  strip 
Jim  t)  <5,  and  g(r)  decays  exponentially  to  zero  as  f Re  t|  -*■  ~ 
uniformly  within  the  strip. 

Similarly,  for  the  same  6,  the  kernel  K(t.t')  is  analytic, 
as  a  function  of  both  t  and  t',  in  the  region 

0  »  |  j  1m  t)  <  8  y  j Im  t*  )  <  6 

In  addition,  K(t,t'}  decays  exponentially  to  zero  as  |Re  t|  ■+•  «, 
uniformly  in  |lm  t|  <5,  provided  t'  is  held  to  a  constant  within  the 
strip  Jim  r ' |  <5.  Similar  behavior  holds  as  a  function  r'  when  r 
is  constant. 


The  properties  described  in  Lemma  1  will  be  used  in  Chapter  V 
in  order  to  develop  an  iterative  solution  to  equations  (4.59-60). 

Once  again  the  author  has  not  been  able  to  derive  a  closed 
form  solution  for  the  unknowns  represented  by  f(t),  except  in  the 
special  cases  n  *  0  or  <p  *  0,  w/2.  However,  the  formulation  of  the 
problem  in  terms  of  Fredholm  integral  equations  of  the  second  kind 
allows  that  vast  wealth  of  knowledge  [37,38,53,54]  regarding  these 
equations  to  come  into  play,  particularly  for  iteratively  generating 
convergent  series  solutions. 


In  the  next  chapter,  an  iterative  series  solution  is  described 
based  upon  the  method  of  successive  approximations  (Neumann  series 
expansion)  from  linear  operator  theory,  and  the  convergence  of  the 
series  is  discussed. 


CHAPTER  V.  THE  METHOD  OF  SUCCESSIVE  APPROXIMATIONS 


5.1  Review  of  Linear  Operator  Theory 

Before  deriving  a  series  solution  of  the  integral  equations 
(4.59-60),  it  is  appropriate  to  review  some  aspects  of  linear 
operator  theory  which  ensure  the  convergence  of  the  series.  Most  of 
the  results  are  taken  directly  from  References  55  through  57. 

A  complex  linear  vector  space  is  a  set  X  together  with  the  set 
of  complex  numbers  Z,  such  that  for  all  x,y€  X  and  a.bG  Z: 

(1)  x+y  *  y  +  x€X 

(2)  ax€X 

(3)  a(x  +  y)  *  ax  +  ay  and  (a  +  b)x  =  ax  +  bx 

(4)  a(bx)  *  (ab)x 

(5)  1  •  x  -  x. 

A  normed  linear  vector  space  is  a  linear  vector  space  X 
together  with  a  function  ||x||  on  X,  such  that  for  all  x,y€  X 
and  a  6  Z: 

(1 )  ]  | x 1 1  >0  and  1 1 x j  |  =  0  if  and  only  if  x  s  0 

(2)  1 1 ax 1 1  =  ( a (  .  ||x|| 

(3)  1 1*  +  y| I  <||*||  ♦  l|y||. 

By  defining  the  function  g(x,y)  a  ||x  -  y||,  for  x,y€  X,  the  normed 
linear  vector  space  X  is  made  a  metric  space  with  metric  g. 
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If  X  Is  a  metric  space  with  metric  g  defined  as  above,  then 
a  Cauchy  sequence  is  a  sequence  xn  in  X  such  that  for  each  real  6  >  0, 
there  exists  a  positive  integer  N  for  which 

g(xn»xm)  ■  llxn  “  <  &  whenever  n.ra  >  N  . 

A  metric  space  X  is  said  to  be  complete  if  for  every  Cauchy  sequence 
xn  in  X  there  exists  a  yG  X  having 

y  *  lim  xn  . 
n  -*■  ® 

A  complete  normed  linear  vector  space  X  is  called  a  Banach 

space. 

If  X  is  a  Banach  space,  then  a  linear  operator  T  is  a  function 
from  X  onto  X  which  Is  linear,  l.e., 

T(ax  +  by)  =  aTx  +  bTy 

for  all  x,y6  X,  a,bG  Z.  The  linear  operator  T  on  the  Banach  space 
X  is  said  to  be  bounded  If  these  exists  some  M  G  [0,®)  such  that 

1 1 Tx J  |  <  M||x||  for  all  xG  X  . 

The  norm  of  the  bounded  linear  operator  T  is  defined  as 


-93- 


xG  A 

where  A  »  {x  €  X|  x  t  0} . 

If  T  is  a  bounded  linear  operator  on  X,  it  is  possible  to  define 
the  nth  iterated  operator  Tn,  n  ^  1,  by  the  inductive  equations 

TJx  =  Tx 
T2x  =  T(Tx) 

Tn+1x  =  T(Tnx) 

for  all  xG  X.  It  is  easily  shown  that 

l|Tn|!  <  | | T I  I "  • 

With  this  background  it  is  possible  to  present  the  following 
theorem. 

Theorem  3: 

If  T  is  a  bounded  linear  operator  on  the  Banach  space  X, 
with  ||T||  <  1,  then  the  series 

ao 

S' 

n=i 

converges  in  operator  norm  to  a  unique  linear  operator  Tq.  Furthermore, 
the  operator 

L  =  I  +  T Q  =  I  +  t  I" 

n=i 
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is  defined  on  X,  where  I  is  the  identity  operator  lx  =  x,  and  from 
| | Tn | |  £  | |T| | n  along  with 


(I  -  T)U  +  t  tP) 
P=1 


=  I  .  T"+1  =  (l  +f  TP)(I  -  T)  , 

p=i 


it  follows  that 

(I  -  T)L  =  I  =  L(I  -  T) 


and  hence  I  -  T  is  one-to-one  onto  X  with  bounded  inverse 

CD 

(I  -  T)"1  =  L  =  I  +  ^  T"  ' 

n=i 

Theorem  3  allows  a  solution  to  the  Fredholm  integral  equations 
(4.59-60)  to  be  formulated  as  a  convergent  series. 

5.2  Series  Solution  to  the  Fredholm  Integral  Equations 

The  general  form  for  the  Fredholm  integral  equations  can  be 
written,  as  before, 

GO 

f(r)  -  A J  K(t,t' ) f ( t ’ )dx'  =  g(x)  ,  (5.1) 

o 

where  f(t)  represents  the  desired  unknown.  Using  operator  notation, 
(5.1 )  can  be  written  as 

(I  -  T)f  =  g  ,  (5.2) 


where  the  linear  operator  T  is  given  by 
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Tf  =  xj  K{T,t')f(x')dr'  .  (5.3) 

o 

It  is  desirable  to  select  a  Banach  space  X  such  that  T  is  a 
bounded  operator  on  this  space.  The  results  of  Theorem  3  can  then 
be  applied  to  equation  (5.2)  to  generate  a  series  solution  for  f(x). 

In  order  to  specify  the  Banach  space  X  it  is  beneficial  to 
recall  the  properties  of  f(x)  and  g(x),  as  they  must  belong  to  the 
space.  In  particular,  from  Chapter  IV  and  Appendix  A, 

1.  f(x),  g(x)  are  analytic  in  a  strip  |lm  x|  <  6,  6  >  0  and 

2.  f(x),  g(x)  +  0  as  | Re  x|  -*■  »  uniformly  within  the  strip 

for  <j>0  -  2ip  >  0,  implying  iji  <  ir/2. 

If  the  region  -°°  <  x  <  ®  is  denoted  by  a,  then  from  (1)  and  (2)  above 

it  follows  that 

sup  |  f(x)  ]  <  00  , 

x  e  a 

and  similarly  for  g(x),  i.e.,  f  and  g  are  bounded  in  a. 

Therefore,  consider  the  linear  vector  space  X  defined  as  the  set 
of  bounded  functions  of  x  analytic  in  the  strip  ( Im  x |  <  6  and 
hence  continuous  on  the  region  a.  The  norm  of  f(x)€  X  is  defined  as 

||f||  =  sup  | f(t) |  .  (5.4) 

t  Ga 

From  a  theorem  of  topology  (see  for  example,  pp.  84,  108,  216  of  [57]), 
it  follows  that  X  is  complete,  and  hence  X  is  a  Banach  space. 
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In  addition,  from  (5.3)  and  the  definition  of  the  norm  it  follows 
that  Tx€  X  for  all  xG  X,  and  thus  T  is  indeed  a  linear  operation 
from  X  into  itself. 

With  the  norm  as  defined  in  (5.4),  it  is  possible  to  show  that 
the  operator  T  is  bounded  for  the  kernels  K(t.t')  occurring  in 
(4.59-60).  In  particular. 


| | Tf I !  =  SUP  |X  f  K(T,T,)f(x')dT,| 

t€E  ol  0 

00 

£  sup  I  X  I  f  |K(T,T‘)|dT'  sup  I f (t ) | 

d  o  t 

00 

1  |x|  1 1  f  1 1  sup  f  |  K(t,t'  ) |  dr '  . 
t€<*  o 

Consider  the  various  kernels  of  (4.59-60). 


(5.5) 


E-Polarization: 


j*  I  Ke(T  ,t  *  )|  dx ' 


ch  ipT 


s 


Ch(lT  -  i|i)t' 
ch  ITT  '  +  ch  ITT 


dt' 


_  ch  IpT  Sh(tT  -  Ip ) T 
sh  iTT  sin  if> 


1 

1 


CSC  if/ 1  1 


sh  ijrr 


ch  4>t 


5 

0 

oc 

l 


_ shU  -  u>)t  1  a  ( 

ch  ITT  '  +  ch  ITT 


ch(*  -  ’to?'  dT. 

ch  ITT'  +  ch  1TTT 


sh ( it  -  2^)t 
sh  ITT 
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It  is  now  possible  to  derive  a  set  of  bounds  for  the  norms  of 
the  various  operators.  Let  the  operator  T  associated  with  the  kernel 


K6(t,t')  be  denoted  by  Tfi,  and  so  on.  From  the  definition  of  the  norm 

of  a  linear  operator,  it  follows  that 

1 1  Te *Tq  1  1  <  \  M  esc  * 

(5.8) 

IIt;.t;h  <  J  |x*|  csc * 

(5.9) 

where 

x  =  -1 

n 

x-  -  -  a 

2  * 

In  order  to  make  use  of  Theorem  3,  it  is  required  that 

and  hence  from  (5.8-9) 

l|T||  <  1, 

In |  >  esc  <|i  ;  for  Te,TQ  , 

(5.10a) 

|n|  <  sin  *  ;  for  T^,T^  . 

(5.10b) 

Equation  (5.10)  provides  sufficient  conditions  on  n  for  the  convergence 
of  the  series  in  Theorem  3. 

Specifically,  the  solutions  to  (4.59-60)  in  the  Banach  space  X 
can  be  written  as  follows: 


V 
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ffi(T)  =  96(t)  +  ^  Tjjge(x) 

n=i 

00 

(5.11a) 

f0(r)  -  g0(T)  +  2  Tjg0CT) 

n=i 

(5.11b) 

where 

Tn  g 
e.o^e.o 

n  » 

■  in)  J  Ke,o(,-,,)3e,o(,')  dT' 

0 

and 

}  n |  >  esc  i>  . 

H-Polarizatlon: 

00 

f;<x)  -  g^tx)  *  ^  T4\<t) 

n=i 

00 

(5.12a) 

fi(,)  ■  «;<'>♦  2 

n=i 

(5.12b) 

where 

1  t' n  g' 

e,oae,o 

n  » 

=  (-r) 

o 


and 

| n |  <  sin  t|i  . 

TT 


The  Iterated  kernel  Kn(T,x')  Is  defined  via 
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KMt.t')  *  K(x  ,t  1 ) 

00 

Kn+1(T,T')  a  j  dT"  ,  n>2  . 

0 

Several  observations  regarding  the  solutions  (5.11-12)  can  be 

made : 

1.  For  E-polarization  the  solution  is  in  the  form  of  a  power 
series  in  n*1 »  and  is  convergent  outside  the  disc  |n|  5  esc  <|>.  For 
H-polarization,  the  series  is  in  powers  of  n,  convergent  within  the 
disc  | n |  <( 2^/tt) sin  ip. 

2.  It  may  be  possible  to  extend  the  regions  of  convergence  of 
(5.11-12)  beyond  the  limits  expressed  in  (1)  above,  since  the  limits 
are  simply  sufficient  conditions  for  the  convergence  of  the  series. 

3.  The  two  regions  of  convergence  for  E-  and  H-polarization 

do  not  overlap,  and  hence  cannot  be  used  simultaneously  to  solve  the 
problem  with  an  arbitrarily  polarized  incident  field. 

4.  The  series  fail  to  converge  for  ~  0,  at  least  In  the 
sufficient  sense.  This  is  consistent  with  Senior's  results  [10] 
which  Indicated  that  for  the  half-plane  an  expansion  in  terms  of 

n  also  contained  contributions  of  order  n  In  n,  and  thus  could  not  be 
expressed  as  a  simple  power  series. 

5.  For  E-polarization,  in  the  limit  n  -*■  ®,  the  unknowns 

ffi,  fQ  -*■  0.  This  is  consistent  with  the  fact  that  for  n  ■*■  »,  the 
wedge  ceases  to  exist,  and  hence  the  scattered  field  is  zero. 


Likewise,  for  H-polarlzatlon  and  n  a  0,  the  result  Is  f' ,f'  =■  gi *9A » 

the  known  solutions  for  the  perfectly  conducting  wedge. 

6.  The  series  (5.11-12)  are  uniformly  convergent  for  t  €  (-»,») 

within  the  bounds  on  n  by  virtue  of  the  norm  (5.4).  When  applying 

the  Inverse  transform  to  the  unknowns  fe,fQ,  etc.,  it  Is  therefore 

possible  to  exchange  the  order  of  Integration  and  summation.  This 

allows  an  approximate  solution  to  be  generated,  wi th  each  successive 

terra  increasing  the  accuracy  of  the  approximation. 

The  final  forms  for  the  unknown  transformed  fields  (u,v)  are 

given  by: 

E-Polarization: 

I  -i(ir/2)v 

Ui(v)  =  - -  [fa(-1v)cos  v(ir  -  4» )  +  f  (-lv)sln  v(ir  -  <J>)  +  21c] 

2v  sin  vir 

(5.13a) 


u2(v) 


e-1(Tr/2)v 
2v  sin  vtt 


COS  v(n  -  t|i) 
COS  VTp 


cos  v<t> 


*  st"  *j>  ;  t)  S1nv«,  ♦  21c  1 

o  sin  vxp  J 

(5.13b) 


H-Polarlzation: 


-  -1(ir/2)v  r  1 

VjM  *  2\>  sin  vtt  [  fe<_1v)cos  “  ♦)  +  fo(-iv)sin  v(tt  -  ♦)  +  21dJ 


*.<»>  *  - [  V-M  ♦>  -  v* 


(5.14a) 


+  f'(-lv)  C0S  Sin  v$+  21  d  1  , 

C  COS  V<p  2  J 


(5.14b) 


where  the  substitution  t  ■  -iv  has  been  used. 
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Recall  that  to  invert  (5.13-14)  directly  requires  restraining 
4  to  regions  where  the  geometrical  optics  fields  are  zero.  Otherwise, 
the  transforms  of  these  fields  must  be  subtracted  prior  to  inversion 
(see  equation  (4.48)  and  Reference  23).  Having  already  assumed  the 
condition  $  -  2ip  =  e  >  0,  the  appropriate  interval  on  $  for  which 
the  geometrical  optics  field  is  zero  is  given  by 

[ ir  ~  4> {  <  c  (5.15) 

in  Region  1,  i.e.,  for  u  ,v  .  In  Region  2,  no  such  interval  exists. 

li  ~ 

Since  it  is  a  difficult  task  to  determine  the  geometrical  optics 
field  in  the  interior  of  the  wedge,  especially  for  small  values  of 
iji,  the  half-angle  of  the  wedge,  one  must  first  determine  (u^.v^) 
from  (ui>vi)  in  the  Interval  (5.15),  and  then  analytically  continue 
the  result  to  the  surface  of  the  wedge.  From  the  boundary  conditions 

it  is  then  possible  to  determine  (u  ,v  )  on  the  surface  of  the  wedge, 

2  2 

and  hence  everywhere  inside  Region  2.  A  more  straightforward 
procedure  is  to  use  Jones'  inversion  formula  (4.4c),  and  to  deform 
the  contour  of  Integration  in  such  a  manner  which  allows  e  to  go  to 
zero.  This  procedure  Is  outlined  in  [22],  In  fact,  Jones'  formula 
can  be  used  to  directly  Invert  all  of  (5.13-14)  without  need  for  a 
restriction  such  as  (5.15). 

It  would  next  be  desirable  to  complement  the  series  solutions 
(5.11-12)  with  similar  results  for  small  n  (E-polarlzation)  or 
large  n  (H-polarizatlon) .  Such  an  approach  Is  considered  in  the 
next  section. 
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5.3  Operator  Theory  for  the  Difference  Equations 

For  the  purposes  of  this  discussion  it  is  sufficient  to 
consider  one  of  the  Integral  equations  (4.59-60),  say  (4.59a). 

It  is  easily  shown  that  the  corresponding  difference 
equation  (4.39a),  which  was  derived  under  the  assumption  that 
Ffi(v)  was  analytic  in  a  strip  | Re  v|  <  5,  where  5  >  1,  can  be 
obtained  from  (4.59a)  under  the  same  assumption.  By  replacing  t 
%t  alternately  t  +  i  and  t  -  i  in  (4.59a),  it  follows  that: 


%«'  ♦  }  J  *£  * 11  •  «.<<  *  <> 


(5.16a) 


ge(t  -  D  . 


(5.16b) 


Equation  (4.59a)  Is  valid  for  | Im  t [  +  |Im  t'|  <1,  which  becomes 
|Ira  r|  <1  since  t'  Is  real.  It  then  follows  that  (5.16a)  is  valid 
for  -2  <  | Im  t |  <  0,  while  (5.16b)  is  valid  for  0  <  | Im  T|  <2. 

From  therelation  (4.58)  between  fe(t)  and  Ffi(v),  where 
v  ■  i-r,  the  requirement  that  F  (v)  Is  analytic  in  | Re  v|  <  5 ,  6  >  1 , 
Implies  f  (t)  Is  analytic  in  |Im  t|  <  6.  This  behavior  allows  one 
to  consider  the  limits  as  |Im  t|  -*■  0+,0'  In  (5.16a,b),  respectively. 
From  the  formulas  of  Plemelj  (see  for  example,  p.  232  of  [58]), 
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Ira  t  o  ,o  o 


-  *  ch  ti/( t  *  i)chU  -  ii>)t  -  /  \  +  ,  ch(it  -  dil-c1  ch  \1>(t  -  i ) 

+  sn  tit  re  v“  /  >  ch  ifT  -  ch  m 


•  fe(t’>  dr*  (5.17) 

Division  of  (5.16a,b)  by  ch  <|»(t  *  i),  taking  the  limit  as  |Im  t| 

0+,0“,  respectively,  and  making  use  of  (5.17)  leads  to 


fe(t  +  i) 
ch  i>( 


-1 _ 1_  -  li  cft(?  Z,#)*  f  (t)  +  ^  T  -*).V - f  (T')dT1 

T  +  i)  n  sh  iTT  Te'T'  3  ch  iit  -  ch  irr  Tev  ' 


96(t  *  i) 
ch  (t  +  i ) 


+  21  f  (t)  +  C  ch(u  -  f  (T-)dT. 

ch  ^(t  -  i )  r)  shit  e  j  ch  itt  -  ch  irx  e 


9e(^  -  *) 
ch  V(t  -  i ) 


Subtracting  the  former  from  the  latter  and  Inserting  the  expression 
for  geCt )  gives 


ffi(t  -  D  ffi(T  +  i) 
ch  iJi(t  -  TT  ch  $ (t  +  i ) 


chU  -  fr) 

sh  ITT 


f.<’> 


16  ch(ir  -  4>0)t  ch  !|»t 
n  sh  ttt 


(5.18) 
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which  easily  reduces  to  (4.39a)  by  letting  t  3  -iv,  fg(-iv)  ■  Fg(\>), 
and  recalling  that  Fg(v)  is  an  even  function. 

The  above  derivation  suggests  an  interesting  result  from 
operator  theory.  In  (5.2)  the  Integral  equation  was  written  as 

(I  -  Tfi)fe  3  ge  ,  (5.19) 

with  Tg  being  an  integral  operator.  Equation  (5.18)  is  then 
equivalent  to 


Tg1  ECI  -  Tfi)fe]  -  (T;1  -  l)fe  -  T;1ge  (5.20) 


where  the  operator  T”1  is  the  left  Inverse  of  T  ,  provided  the  space 

C  C 

of  functions  Is  restricted  to  those  which  are  analytic  in  a  strip 

| Im  t|  <  6,  6  >  1.  It  is  also  possible  to  show  that  T"1  is  a  right 

Inverse  under  the  same  restrictions.  Explicitly,  T!1  Is  defined  by 

0 


T;if. 


f,(T  -  1)  f,<T  ♦  D 

ch  tp  ( x  -  i)  ch  ip  ( x  +  i ) 


(5.21) 


where  f  Is  even  and  analytic  in  | Im  t|  <  6 ,  6  >  1 .  Obviously, 
Tg1  is  a  difference  operator. 

Equation  (5.20)  can  be  written  as 


u  -  ?>f. 


•T-1g 
e  ae 


(5.22) 
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whlch  leads  one  to  ask  If  the  results  of  Theorem  3  can  be  used  to 

generate  a  series  solution  for  f  In  terras  of  the  Iterated  operator 

(T"1 )n,  which  Is  essentially  a  series  In  powers  of  n.  Such  a  series 

would  Indeed  complement  the  earlier  result,  equation  (5.11a),  which 

-l 

expresses  ffi  as  a  power  series  In  n  • 

However,  the  author  has  not  been  able  to  develop  a  suitable 

bound  for  the  norm  of  T'1,  and  hence  Is  not  able  to  take  advantage 

of  the  results  of  the  theorem.  Nonetheless,  with  a  suitable  choice 

of  a  Banach  space  of  functions  and  operator  norm.  It  may  be  possible 

to  prove  the  convergence  of  a  series  in  T”1.  Indeed,  the  author 

has  shown  that  the  first  two  terms  of  such  a  series  produce  the 

correct  geometrical  optics  fields.  Furthermore,  when  n  ■  0,  the 

series  reduces  to  a  single  term  which  is  the  correct  solution  for  the 

perfectly  conducting  wedge.  Thus  It  appears  that  the  series,  even  if 

It  Is  not  convergent,  may  provide  an  asymptotic  representation  as 

n  -*■  0.  Similar  conclusions  regarding  the  other  operators  T  ,  T', 

T'  can  be  made, 
o 


CHAPTER  VI.  COMMENTS  AND  CONCLUSIONS 


The  electromagnetic  scattering  of  an  arbitrarily  polarized  plane 
wave  normally  incident  upon  a  resistive  wedge  has  been  formulated  via 
a  pair  of  related  function-theoretic  techniques,  the  method  of 
Maliuzhinets  [12-14]  and  the  Kontorovich-Lebedev  transform  method 
[20-22],  both  of  which  have  been  successfully  used  to  solve  scattering 
problems  in  a  single  wedge-shaped  region. 

The  goal  of  the  author  was  twofold:  obviously  to  find  an  exact 
solution  to  the  resistive  wedge  problem,  and  secondly,  to  point  out  the 
similarities  between  the  two  function- theoretic  techniques  as  well  as 
the  complexity  that  arises  in  applying  them  to  a  two-region  problem. 

With  regard  to  the  latter  of  these  goals,  the  author  has  shown 
that,  as  with  past  applications,  both  methods  lead  to  a  set  of  difference 
(or  functional)  equations  for  the  various  unknowns.  However,  If  these 
methods  are  simply  applied  directly,  as  in  the  single  region  problem, 
the  presence  of  nonzero  fields  in  two  regions  produces  a  coupling  of  the 
unknowns  in  the  functional  equations  (see  equations  (3.9-10)  and 
(4.31-32),  for  example).  By  means  of  a  decomposition  into  symmetric 
and  antisymmetric  components,  the  equations  were  successfully  uncoupled, 
although  they  remained  sufficiently  complex  to  prevent  the  determin¬ 
ation  of  closed  form  exact  solutions. 

One  reason  for  this  shortcoming  is  the  lack  of  a  systematic 
technique  for  solving  difference  equations  of  order  greater  than  one. 
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To  this  end,  the  author  has  developed  a  novel  procedure  which, 
under  less  restrictive  conditions,  replaces  the  difference  equations 
of  the  K-L  method  with  Fredholm  integral  equations  of  the  second 
kind.  This  allows  a  large  number  of  techniques  to  come  into  play 
for  the  determination  of  exact  and  approximate  solutions. 

One  such  technique,  the  method  of  successive  approximations, 
is  shown  to  lead  to  uniformly  convergent  power  series  solutions  for 
certain  values  of  the  normalized  resistivity  n.  In  particular,  for 
E-polarization,  the  series  converge  when  n  is  large,  while  for 
H-polarization,  they  converge  when  n  is  small.  This  behavior 
prevents  the  use  of  both  polarizations  simultaneously  for  generating 
the  solution  to  an  arbitrarily  polarized  incident  field. 

Although  an  attempt  to  alleviate  this  problem  was  made  by 
demonstrating  that  the  integral  operator  and  difference  operator 
were  inverses  of  each  other  under  restricted  conditions,  the  author 
was  not  able  to  bound  the  difference  operator  over  the  space  of 
functions  being  considered,  and  hence  could  not  prove  a  series 
generated  with  this  operator  would  converge.  Nonetheless,  the  series 
does  exhibit  proper  behavior  under  certain  circumstances,  leaving  open 
the  question  as  to  whether  a  bound  for  the  operator  does  exist. 
Certainly,  this  is  an  area  open  to  future  work.  A  parallel  effort 
Investigating  more  accurate  bounds  for  the  integral  operators, 
thereby  extending  the  radius  of  convergence  of  the  series,  is  also 
worth  pursuing. 

In  addition,  the  author  is  hoping  to  communicate  the  results 
of  this  work,  particularly  the  functional  equations  (3.32-33)  and 


(4.39-40),  to  several  authorities  on  functional  equations  for  their 
suggestions  regarding  possible  solutions. 

Finally,  the  reader  may  find  this  work  conspicuously  lacking 
in  numerical  or  graphical  results.  This  is  not  an  accident,  as  the 
purpose  of  the  author  was  to  explore  the  possibility  of  extending 
proven  function-theoretic  techniques  toward  finding  an  exact  solution 
for  the  resistive  wedge,  and  in  the  process,  to  illustrate  the  nature 
of  the  complications  that  arise  when  these  techniques  are  applied 
to  a  two-region  problem.  The  author  believes  that  the  presentation 
of  numerical  data  based  on  an  approximate  solution  is  not  consistent 
with  these  goals,  and  hence  has  reserved  this  area  for  pursuit  in 


the  future. 


APPENDIX  .  THEOREMS  FOR  THE  METHODS  OF  MALIUZHINETS  AND  THE 

KONTOROYICH-LEBEDEV  TRANSFORM 

A.l  The  Method  of  Maliuzhinets 

The  material  in  this  section  is  derived  entirely  from  the 
works  of  Mai iuzhinets  [13-15,51],  The  basis  for  his  method  is 
the  representation  of  a  function  S(p,<{>)  of  the  polar  coordinates 
(p,<j>)  in  the  form  of  a  Sommerfeld  integral,  viz. 

S(p,<f>)  =  2rf  J  elkpC0S<x  s(a  -  <f> )dot  (A.l) 

Y 

where  y  is  a  contour  in  the  complex  a  plane  consisting  of  two  loops, 

Y  and  y2.  The  contour  approaches  infinity  within  the  shaded  regions 
where  the  real  part  of  ikpcosa  is  negative  for  positive  real  k,  as 
shown  in  Fig.  3.1,  and  reproduced  for  convenience  in  Fig.  A.l. 

In  [14],  Maliuzhinets  establishes  conditions  for  the  existence 
of  a  unique  solution  s(a)  to  the  integral  equation  (A.l)  for  S(p,4>) 
satisfying  certain  boundedness  conditions. 

Theorem  A.l : 

Let  M,a,b,c,d  be  positive  numbers,  and  let  e,m  be  numbers 
satisfying 

0  <  e  <  it 

[arg  m|  <  | 
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Let  F(p)  be  a  function  satisfying  the  inequality 

I  F(p )  |  <  M[Pr1+a  eb*pl 

for  positive  values  of  p,  and  also  in  the  entire  region 

0  <  |  p  |  <  « 

j  arg  p  |  <  ,  0<e1<ir  , 

where  this  function  is  analytic.  Consider  the  integral  equation 

F(p)  -  2“T  J  emp,cos“  f (a)  da  ,  (A. 2) 

y' 

where  the  contour  y'  is  made  up  of  two  loops,  yj  and  y^  .  The  loop 

y  consists  of  the  two  half  lines 
1 

Re  a  =  arg  m  t  j 

Im  a  >_  d 

and  the  line  segment  Im  a  =  d.  The  loop  y'  is  symmetric  with 

2 

respect  to  y^  about  a  =  0  (see  Fig.  A. 2). 

Then  there  exists  one  and  only  one  solution  f(a)  to  (A. 2) 
which  is  analytic  on  and  within  the  contour  y'  except  at  infinitely 
distant  points,  and  which  satisfies  the  additional  constraints 
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(1-a  ) | Ira  a | 

lf(a)|  <  M  e  1  ;  M  ,a  >  0 

2  12 

f(ot)  =  -f(-a)  . 


This  function  is  represented  by  the  integral 


f(«)  - 


m  si 
2 


—  I 


F(p)e 


-mpcosa 


da 


(A. 3) 


(A. 4) 


for  Re(m  cos  a)  >  0.  For  this  function,  a  =  a. 

l 


Proof: 


In  view  of  the  fact  that  f(a)  is  odd,  (A. 2)  can  be  written  as 


F(p)  =  V  f  empC0Sa  f(o)  da  (A. 5) 

in  j 

By  making  a  change  of  variables,  W  *  en^ar®  n^cosa,  and  defining 


g(W)  -  .  24^_e-1(ar3  "> 
3  sin  a 


(A. 6) 


equation  (A. 5)  becomes 

F(p)  *  2^T  J  el™  lpW  g(W)  dW  ,  (A. 7) 

r 

where  r  is  the  image  of  yj  under  the  transformation  of  variables.  The 
contour  r  Intersects  the  real  axis  between  zero  and  ch  d,  and  coincides 
at  infinity  with  the  rays  arg  W  *  ±(e  +  j  ),  as  shown  in  Fig.  (A. 3). 
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The  function  g(W)  is  analytic  to  the  right  of  r,  which  coincides 
with  the  interior  of  the  loop  where  f(a)/sin  a  is  analytic. 

From  (A. 3)  it  follows  that 


|g(W)|  <  4MjW|“ai  ,  |W|  -  =0  ,  |argW|<£+!  .  (A. 8) 


Multiplying  (A. 7)  by  e”^pW,  where  Re  W  >  ch  d,  and  integrating 
over  (0,«>)  on  p  gives 

00  00 

^  F(p)e‘lmlpW  dp  =  ~  j  dp  J  elmlp(WrW)  g(Wi)  dWi  ,  (A. 9) 

0  or 

where  the  integration  is  taken  as  a  limit  at  the  endpoints.  By  virtue 
of  (A. 8)  and  the  fact  that  ReCV^-W)  <  0,  the  order  of  integration  may 
be  exchanged,  and  the  limit  taken,  which  leads  to 


J  F(P)e'MpW  dp 


_  i  ,  f  SOU 

2wTfmi  J  vnnr  dwi 

r  1 


(A. 10) 


where,  thanks  to  (A. 8),  the  integral  in  (A. 10)  has  been  evaluated  from 
the  theory  of  residues.  Transforming  W  back  to  a, 

00 

S  F(p)  dp  -  -  MfL 

0 

which  is  the  desired  result  given  in  (A. 4). 

In  order  to  show  that  the  assumptions  regarding  the  analyticity 
of  f (a) ,  along  with  the  conditions  (A. 3),  are  indeed  true,  consider 
equation  (A. 10)  as  a  definition  of  the  function  g(W). 


V 
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Because  of  the  conditions  that  F(p)  be  analytic  for 
|arg  pj  <  e  ,  it  follows  that  the  integration  for  p  >  c  in  (A. 10) 
may  be  displaced  to  an  arbitrary  half  line  |arg  p|  =  e  >  e  (by  choice), 
from  which  it  follows  that  g(W)  is  analytic  for  large  jWj,  provided 
|arg  Wj  <  e  +  n/2.  Defining  f(a)  via  (A. 6),  one  then  concludes  that 
for  sufficiently  large  d,  f(a)  is  analytic  in  Im  a  >_d,  | Re  a  -  arg  m| 

<  e  +  ir/2 ,  which  corresponds  to  the  interior  of  the  loop  y ' . 

In  addition,  because  it  has  been  assumed  that 

|F(p)|  <  H|B|-Ua  ebM  , 

It  can  be  shown  via  (A. 10)  that 

Ifljp  <  j  | f(p )|  «>  dp 

0 

00 

<  M|mW|’a  j  x  -1+a  e-x[cos(argW)-b/|<nW|]  dx  _ 

o 

Since  the  Integral  is  bounded  for  sufficiently  large  |W|, 

I g(W) I  <  M  |W|"a  ,  M  >  0  , 

2  2 

confirming  (A. 8),  and  since 

= 

sin  a 


-2Wg(W) 
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1t  follows  that 


|  f (a) j  <  M1  e(1-al ) I Intal 


with  ax  =  a,  as  assumed.  Thus  the  proof  is  complete. 
By  means  of  the  substitutions 


F(p)  *  S(p  ,4) ) 


f(a')  =  y  [s(a*  +  tt  -  4» )  -  s(-a'  +  it  -  4> ) 3 


M  =  ik  (implying  Im  k  >  0) 


and  replacing  a'  «  a  -  ir,  equation  (A. 5)  becomes 


S(p,<f>)  = 


J 


kpcosa 


[s(a  -((>)-  s(-a  +  2ir  -  <(>)]  do 


From  Fig.  A.l  and  the  relation  between  the  contours  y  and  y  , 

1  2 


the  expression 


S(p,$)  s  2rf  J  eikpCOSot  s(o  -  $)  do 


is  obtained,  which  agrees  with  (A.l). 

An  important  result  of  Theorem  A.l  is  that  an  odd  solution 
to  the  equation 


jU-  f  fl.')e”<,C0Sa'  da'  .  0 
y' 

must  vanish  identically;  consequently,  a  solution  to  the  equation 

j  eikpcosa  s^a  _  ^  .  o  (A. 12) 

Y 

roust  satisfy  the  condition 

s(a  -  $)  =  s(2ir  -  a  -  <t>)  ,  (A. 13) 

i.e.,  the  coefficient  of  the  exponential  in  the  integral  must  be  even 
about  a  »  it.  This  result  is  used  extensively  in  formulating  the 
functional  equations  for  the  unknowns  in  Chapter  III. 

The  justification  of  a  representation  such  as  (A-l)  for  the 
unknown  fields  satisfying  the  Helmholtz  equation 

(v2  +  k2)S(p,d>)  =  0 

Is  discussed  in  considerable  detail  In  [51],  and  will  not  be  reproduced 
here  for  the  sake  of  brevity.  However,  it  is  worthwhile  noting  that  the 
bounds  placed  upon  the  function  F(p )  in  Theorem  A.l  include  the  functions 
S(_p )  satisfying  the  edge  and  radiation  conditions  discussed  In  Chapter 
II.  Furthermore,  from  [51],  the  value  of  S(p,$)  at  p  *  0,  provided 
S(p4)  Is  bounded,  is  given  by 
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S(0  ,<j> )  =  2is(i»)  =  -2is(-i»)  ,  (A. 14) 

and  does  not  depend  on  as  expected.  The  boundedness  at  p  =  0  Implies 
a  ^  1 ,  which  also  implies,  from  (A. 3),  that  s(i«)  <  »,  consistent  with 
(A. 14) . 

The  representation  (A-l)  is  simply  a  superposition  of  elementary 
plane  wave  solutions  to  the  Helmholtz  equation.  This  is  more  readily 
seen  by  replacing  a  by  a  -  ir  +  <fi,  whereby 

5(0.4.)  •  2JT  ^  e-lkocos(«+*)  s(o  .  „)  (A.I5) 

T* 

where  y,  is  simply  the  contour  y  displaced  an  amount  to  the  right. 

As  discussed  in  [51],  (A. 15)  represents  a  set  of  plane  waves  incident 
from  a  direction  2-n  -  a,  where  a  is  complex,  allowing  for  evanescent 
(decaying)  waves. 

Recalling  the  geometry  for  the  scattering  of  a  plane  wave  by  a 
resistive  wedge,  the  incident  field  is  given  by 

(uV)  =  e-1kpcos(ij>0-<j>)  _  (A. 16) 

Maliuzhinets  shows  in  [51]  that  discrete  plane  waves,  given  by  the 
geometrical  optics  fields,  correspond  to  poles  of  s(a-ir)  lying  within 
the  strip  ir-$<Rea<3ir-$.  From  (A.  16),  it  follows  that  s(a-ir) 
has  a  pole  at  a  3  2*  -  $  with  residue  unity.  Furthermore,  the  field 
(A. 16)  is  the  only  geometrical  optics  field  incident  from  within 
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Region  1  (<|>  <.  $  <.  2*  -  i(») ,  as  all  other  fields  appear  to  emanate 
from  within  the  wedge  (i.e..  Region  2). 

Therefore,  it  may  be  stated  that  s^o)  (representing  the  field 
in  Region  1)  satisfies  an  additional  constraint  for  the  wedge  problem, 
that  is 

(a  -  it)  -  (a  -  2ir  +  <j>Q)  1  (A. 17) 

is  analytic  for  ip  <_  Re  a  <_  2tt  -  .  This  condition  is  stated  in 

Chapter  III  as  a  means  of  uniquely  determining  a  solution  to  the 
functional  equations  derived  for  the  unknowns. 

A. 2  The  Kontorovich-Lebedev  Transform 

The  requirements  for  the  existence  of  the  K-L  transform  and  its 
inverse  are  rigorously  described  in  [21,22,59],  each  of  which  prescribes 
a  set  of  conditions  for  a  function  f(p)  or  f(v)  in  order  for  the 
integral  formulas  to  converge.  The  results  of  [21]  have  been 
reproduced  without  proof  in  Theorem  2,  Chapter  IV  of  this  work. 

Rather  than  repeat  any  of  these  derivations  here,  the  author 
will  simply  state  the  necessary  results,  and  then  apply  them  to  the 
particular  functions  being  considered  in  Chapter  IV  in  order  to 
justify  the  assumptions  made  therein. 

Specifically,  a  set  of  total  fields  (U,V),  representing  the 
solution  to  the  scattering  by  a  resistive  wedge  of  an  E-  or  H-polarized 
Incident  plane  wave,  were  shown  in  Chapter  II  to  have  certain 
boundedness  properties  as  p  approaches  zero  and  infinity. 
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Near  the  vertex  of  the  wedge 

(U,V)  =  0(p5)  ,  6>0  ,  P  +  0  ,  (A. 18) 

which  allows  (U,V)  to  be  at  most  a  nonzero  constant  at  p  =  0.  Further¬ 
more,  by  writing, 


(U,V)  =  (U9  +  Ud,  V9  +  Vd) 


(A. 19) 


i.e.,  as  a  sum  of  geometrical  optics  and  diffracted  fields,  the 
radiation  condition  requires 


(A. 20) 


which  implies  that  (Ud,Vd)  decay  at  least  as  rapidly  as  p  as 
p  •>  In  addition,  several  of  the  inversion  theorems  assume 
Im  k  >  0,  in  which  case,  writing 


(U,V)  =  (U1  +  Us,  V1  +  Vs)  , 


the  scattered  fields  (IIs, Vs)  behave  as  e’bp  as  p  +  »,  where  b  >  0. 
It  is  convenient  to  summarize  these  results  into  a  set  of 


assumptions  regarding  (U , V) : 


1.  There  exist  A,  pQ  >  0  such  that 


|u,y(  <  Apa  ,  6  >  0 


(A. 21) 


for  p  <  p. 


2.  There  exist  B,C  >  0  such  that 


|Ud,Vd|  <  Bp" 1/2  eikp 
|  U9 ,V9 1  <  c 


(A. 22a) 
(A. 22b) 


for  p  >  pQ. 


■I 


Since  the  value  of  |U  ,V  |  is  unity  for  all  p,  it  follows 


that  (US,VS)  satisfy  (A. 21)  above.  Inasmuch  as  the  scattered  fields 


contain  geometrical  optics  terms,  (IIs, Vs)  also  satisfy  (A. 22b) 
In  Chapter  IV,  the  modified  scattered  fields 


Js  +  c  e1kp 


Vs  +  d  e1kp 


where  c  =  -Us  (p=  0),  d  =  -Vs  (p  =0),  were  defined.  This  implies 
that  (u,v)  satisfy  (A. 21)  with  6  strictly  greater  than  zero.  In  a 
manner  analogous  to  (A. 19),  it  is  possible  to  separate  the  geometrical 
optics  terras  from  (u,v),  (as  was  done  in  [23]),  viz 


u  *  u9  +  ud  =  (Usg  +  cg  e1kp)  +  (Usd  +  cd  eikp)  (A. 24) 

and  similarly  for  v.  Note  that  ug,vg  both  satisfy  (A. 21)  with  &  >  0, 
in  addition  to  (A. 22b). 
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Consider  the  K-L  transform  of  (u,v),  given  by 

oo 

(u,v)  =  J  (u.vjH^Ckp)  .  (A. 25) 

o 

Because  H^(kp)  =  0(p"^ev^},  it  follows  from  (A. 21)  and  (A. 22b)  that 
the  integral  (A. 25)  is  uniformly  convergent  for  [ Re  v |  <6,  and  thus 
(u,v)  are  analytic  functions  in  the  strip  | Re  v|  <  <S.  In  addition, 
from  [23]  it  can  be  shown  that  (u^,v^)  approach  zero  as  |  Im  v|  «, 

| Re  v |  <  <5.  Unfortunately,  the  geometrical  optics  transformed  fields 
do  not  exhibit  such  behavior. 

In  particular,  a  general  form  for  (u9,v9)  is 

N 

(u9,v9)  =  Y  an[e'ikpC0Sen  -  eikp]  (4.26) 

n=l 


where  en  =  <}>  -  4>n,  <j>n  is  the  negative  direction  of  propagation  of  the 
plane  wave,  and  an,N  are  constants.  The  transform  is  then  given  by 


N 

(Ug,v9)  =  Za(  (e-iXpcose  .  e^P)H(1  )(kp)  dp 

n=1  n  J  V  0 


n=l 

N 


.  ^  2ia  e"l(ir/2)  1  ~  cos  v(ir  *  gn> 

n=l 


v  sin  vir 


(A. 27) 


Note  that  the  second  order  pole  at  v  =  0  is  cancelled  by  the  second 

order  zero  of  [1  -  cosv^-l^)].  Therefore  (u9,v9)  are  analytic  in 

| Re  v J  <  1.  The  minus  sign  corresponds  to  0  <_  en  <.  2*,  the  plus  sign 

to  -2ir  <  b  <0-  It  is  easy  to  show  that  (u9,v9)  become  unbounded 
-  n  — 

for  ir/2  <  |en|  <  3ir/2  when  v  -1®.  It  is  for  this  reason  that  one 
cannot  perform  the  inversion 


'V 
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•jco 

(u,v)  =  -  2  ^  (u,v)v  JvCkp)  dp 

-i® 


directly.  Instead  it  is  the  integral 

i“ 

(u  -  u9,  v  -  v9)  =  -  ^  ^  (u  -  u9,  v  -  v9)vJv(kp)  dp  (A. 28) 

-i» 

which  must  be  considered.  Jones  avoided  this  problem  in  [22]  by 
defining  the  inversion  via 

i°° 

(u,v)  a  -  lira  Y  ^  eev2  (u,v)vJv(kp)  dp  ,  (A. 29) 

n+  -  i  oo 

e-*-0 

allowing  the  inclusion  of  the  geometrical  optics  terms. 

It  should  be  noted  that  should  there  exist  intervals  in  <j>  for 
which  (u3,v9)  are  zero,  then  the  inversion  may  be  carried  out  for  these 
values  of  <t>,  and  the  result  analytically  continued  for  other  values. 

This  method  neatly  generates  the  geometrical  optics  terms. 

In  Section  4.2.2,  the  derivation  of  the  integral  equations  for 
the  various  unknowns  made  use  of  a  formula  of  the  form 


j  u9H^  }(kp)  dp 
o 


L.  B-U*/2)v 
2k  e 


sin  uir  u9 

COS  hit  +  COS  Vir 


dp 


.(A. 30) 


It  is  desirable  to  explicitly  show  this  result. 


By  inserting  (A. 26)  into  the  LHS  of  (A. 30)  one  obtains 


\  u9Hv(l)(kp)  dp 
o 


V  2an  e ■^ir/'2^v  s1n  vCw  “  I Sn  I ) 

2j  T~  sin  vtt  v  +  sin|e  I 


.  (A. 31) 
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This  result  was  derived  In  Chapter  IV  In  a  slightly  different  manner, 
based  on  the  uniform  convergence  of  the  transform  representation  of 
u*.  The  results  are  equivalent. 
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